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PREFACE 


In  a summer  1976  workshop  at  the  University  of  Hawaii  an  international  group  met 
under  the  auspices  of  the  Deep  Underwater  Muon  and  Neutrino  Detection  Committee  (DUR- 
MAND)  to  discuss  the  detection  of  cosmic  particles  and  their  interaction  byproducts  by  listen- 
ing through  hydrophones  to  the  noises  they  make  in  the  ocean.  A full  proceedings  of  their 
work  was  published  in  early  1977  (cited  as  Ref.  1 in  this  report).  Of  particular  interest  were  the 
conclusions  of  the  Acoustics  Panel,  which  examined  the  acoustic  aspects  of  DUMAND  and 
developed  several  acoustic  models  to  serve  as  prediction  tools.  Although  the  models  are  accessi- 
ble to  advanced  students  of  acoustics  theory,  they  were  not  given  in  sufficient  detail  to  allow 
appreciation  of  their  content.  The  purpose  of  this  survey  report  is  to  place  the  models  in  per- 
spective by  showing  their  origin,  derivation,  and  limits  of  validity. 

The  conclusions  of  the  DUMAND  Acoustics  Panel  covered  many  features  of  the  acoustic 
problem.  Points  most  pertinent  to  this  survey  are  briefly  as  follows: 

• The  cosmic  particle  deposits  energy  in  the  ocean  as  it  decelerates.  The  spatial  distri- 
bution of  this  energy  along  the  path  of  the  resultant  conical  cascade  is  roughly  Gaussian.  The 
length  of  track  varies  from  1 to  10  m,  depending  on  the  energy  fp  of  the  incoming  particle 
and  on  the  nature  of  the  cascade.  For  hadronic  (massive-particle)  plus  electromagnetic 
showers  the  full  width  at  half  maximum  (FWHM)  of  the  Gaussian  distribution  in  the  direction 
of  the  shower  is  about  3.8  m at  lO’^cFto  6.4  m at  10 eK  The  lateral  distribution  (across  the 
shower)  is  not  specified  but  is  expected  to  be  similar  to  purely  electromagnetic  cascades.  For 
purely  electromagnetic  cascades  the  FWHM  is  3.2  m at  lO'^eFto  4.5  m at  10 In  the 
cross  section  the  fraction  (f^)  of  total  energy  deposition  between  r(r  = 0 is  the  axis  of  the 
cascade)  and  (Moliere  length,  which  equals  10  cm  in  deep  sea)  is  = 2.5  r/rj  if  r < < rj. 

• The  ambient  noise  in  the  ocean  is  a minimum  at  25  kHz,  making  it  a good  cnoice  for 
a listening  frequency.  At  25  kHz  the  ambient  noise  power  (in  SI  units)  is  proportional  to  (10 

in  a 1-Hz  band.  Since  a bandwidth  of  some  10  kHz  is  needed  to  resolve  the  expect- 
ed transients  due  to  the  acoustic  pulses  of  the  cosmic  particles,  the  noise  power  in  the  receiver 
is  at  best  proportional  to  (10/z A//n^) ^10“*,  so  that  the  noise  pressure  spectrum  level  is 

V (lO^xA/m) ^10^  =10“^(A/m^)  = 10  (dyne/cm ^ ).  This  is  taken  as  the  lowest 
DUMAND  signal  that  can  be  detected  in  the  band  15  to  25  kHz  say. 

• The  mass  of  ocean  water  required  to  detect  neutrinos  must  be  very  large  because  of 
the  extremely  small  scattering  cross  section  of  neutrinos.  It  is  estimated  that  a volume  of  10 
by  10  by  1 km  ^ is  needed  to  detect  100  neutrino  related  events  per  day,  each  event  having  an 
energy  averaged  over  many  of  10 eV.  The  number  of  hydrophones  needed  to  make  this 
detection  over  the  prescribed  volume  is  estimated  at  10  ^ 

One  gathers  from  the  DUMAND  report  that  further  progress  in  the  prediction  of  noise 
pulses  from  the  deceleration  of  high-energy  particles  in  the  ocean  must  await  more  detailed  ex- 
perimental data,  specially  designed  to  display  the  impulse  nature  of  the  physical  event  and  give 
reasonable  estimates  of  time  duration  and  magnitude  of  heat  generation.  Although  it  has  been 
verified  in  the  laboratory  that  high  energy  particles  do  make  noise  when  they  are  decelerated 
in  liquids,  the  lack  of  good  experimental  data  precludes  any  decision  on  which  model  is  most 
nearly  correct. 
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THEORETICAL  MODELS  AND  NUMERICAL  ESTIMATES 
OF  ACOUSTIC  SIGNALS 

OF  HIGH-ENERGY  COSMIC  PARTICLES  IN  THE  OCEAN 


EXECUTIVE  SUMMARY 

High-energy  physicists  in  the  field  of  cosmic  rays  are  faced  with  a difficult  problem  that  is 
blocking  advance.  They  wish  to  clarify  their  understanding  of  the  universe  by  probing  into  the 
origin  of  neutrinos  which  reach  the  earth  from  outer  space.  In  particular  they  wish  to  con- 
struct a neutrino  telescope  that  will  pinpoint  the  angular  distribution  of  these  particles  in  the 
heavens,  and  chart  their  temporal  and  spatial  fluctuations.  To  date  they  have  been  unable  to 
do  this  in  any  practical  way,  because  the  scattering  cross  sections  of  neutrinos  are  extremely 
small,  even  at  high  energies.  Groups  of  scientists  interested  in  this  problem  have  assembled 
themselves  into  a committee  with  the  name  DUMAND  (Deep  Underwater  Muon  and  Neutri- 
no Detection)  and  have  met  in  the  last  few  years  at  annual  intervals  to  discuss  advances  to- 
ward solution  of  their  problem.  The  most  recent  meeting  was  at  the  U.  of  Hawaii,  in  the  Sum- 
mer of  1976.  They  conclude  the  following  in  their  report:  immense  masses  of  matter  will  be 
needed  to  detect  these  muons  and  neutrinos  at  some  usable  (practical)  rate,  say  a few  tens  or 
hundreds  of  events  per  day.  This  mass  must  essentially  be  free  to  all  at  no  cost.  Clearly  the 
only  practical  resort  to  get  great  mass  is  the  ocean.  Hence  they  visualize  the  neutrino  telescope 
to  be  an  underwater  detector.  Furthermore,  to  avc'd  too  many  non-muon  and  non-neutrino 
events  per  unit  time,  which  could  cause  saturation  in  the  detector  they  conceive  the  underwa- 
ter detector  to  be  deep  in  the  ocean.  Thus  an  immense  structure  of  hydrophones  and  photo- 
detectors estimated  to  be  some  5 to  10  km  square  and  1 kilometer  high  is  proposed,  eventually 
at  a depth  of  5 kilometers,  in  some  favorable  ocean  far  away  from  man-made  activities  which 
create  underwater  noise. 

The  theory  of  the  neutrino  telescope  is  this:  High-energy  particles  entering  seawater  give 
birth  to  cascades  of  ionized  collision  products.  These  cascades  can  be  detected  by  two  effects: 
they  create  a bow  wave  of  light  (Cerenkov  effect)  (similar  to  a surface  "bow  wave"  caused  by  a 
ship)  which  is  visible  as  a flash  of  light,  and  they  generate  an  acoustic  shock  pulse,  theoretical- 
ly detectable  with  hydrophones.  The  source  of  the  energy  in  the  acoustic  effect  is  thought  to 
be  ihe  ionization  energy  loss  caused  by  charged  particles  colliding  with  water  molecules.  The 
energy  ultimately  appears  as  nearly  instantaneous  local  heating  of  the  water,  in  a time  of  about 

1 ps  to  0.01  ns,  over  a volume  of  water  at  the  largest  of  radius  10  nm  but  generally  of  the  di-  ^ 

mensions  of  a few  molecules.  The  efficiency  of  conversion  of  kinetic  energy  of  cosmic  particle  ij 

into  acoustic  energy  is  estimated  to  be  10  or  less.  i 

i 

From  the  point  of  view  of  the  acoustician  the  problem  is  posed  as  follows.  Coming  in  | 

from  the  atmosphere,  and  ultimately  from  outer  space,  is  a continuous  shower  of  high-energy  i 

particles,  charged  and  uncharged,  which  collide  with  molecules  of  air  to  generate  cascades  of  I 
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collision  products,  ionized  molecules,  etc.  In  particular  a large  flux  of  massless,  chargeless  neu- 
trinos enter  the  seas,  where  they  in  turn  collide  with  water  molecules  and  create  additional  cas- 
cades. At  the  instant  these  cascades  are  formed,  they  radiate  transient  pulses  of  sound.  Thus 
the  acoustician  sees  a volume  of  space  randomly  illuminated  with  bursts  of  sound  (firefly 
effect).  The  problem  posed  then  is  to  create  a mathematical  model  of  this  effect. 

Several  models  of  single  cascades  were  initially  made  under  the  auspices  of  DUMAND. 
These  can  be  grouped  into  thermoelastic  models  and  bubble  models.  In  the  group  of  thermoe- 
lastic models  are  a heated-rod  model,  a heated-spot  model,  and  a heated-filament  model.  In  the 
heated-rod  model  the  cascade  is  considered  to  be  a cylinder  of  water  typically  3 to  10  m long 
and  1 to  12  cm  in  diameter,  heated  in  a time  of  1 /zs  to  1 ns,  and  coherently  radiating  a tran- 
sient whose  spectrum  is  roughly  from  25  kHz  to  several  hundred  kilohertz.  The  radiation  pat- 
tern is  thought  to  be  pancake  shaped,  that  is,  sharply  narrow  broadside  to  the  axis  of  the 
cylinder.  In  the  heated  spot  model  a spherical  volume  of  water  of  radius  10  nm  or  smaller  is 
heated  in  0.1  ns  to  0.01  ps.,  radiating  a spherical  shock  wave  followed  by  a tail  thought  to  be 
due  to  a heat  diffusion  effect.  In  the  heated-filament  model  the  local  heating  is  confined  to  a 
large  number  of  fine  filaments  approximately  1 ftm  in  diam.  Each  filament  radiates  a separate 
acoustic  wave,  but  all  the  filaments  are  considered  to  radiate  coherently.  In  the  group  of  bub- 
ble models  the  principle  effect  is  modeled  as  local  boiling  of  the  fluid  medium  into  bubbles, 
with  subsequent  collapse  of  these  bubbles  giving  rise  to  acoustic  transients.  Not  too  different 
are  models  which  treat  the  phenomenon  as  local  explosions,  as  if  made  of  TNT,  the  equivalent 
weight  being  0.01  fig  of  TNT  or  less. 

Both  the  thermoelastic  models  and  the  bubble  models  have  been  widely  discussed.  They 
are  reviewed  in  detail  in  the  body  of  this  report,  and  a number  of  calculations  have  been  made 
to  illustrate  the  character  of  the  predictions  that  can  be  made  with  their  use. 

Experiments  conducted  to  date  in  laboratory  mockups  [1]  show  that  the  transient  acoustic 
pulse  due  to  high-energy  particles  partakes  of  the  nature  of  a wave  called  a bipolar,  with  an  ini- 
tial compression  wave  followed  by  a final  rarefraction,  as  shown  in  the  following  sketch: 
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This  type  of  acoustic  signal  is  predicted  by  a generalized  thermoelastic  model  developed  in  the 
body  of  this  report.  The  model  is  obtained  by  solving  the  linear  acoustic  wave  equation  driven 
by  a heat  source,  namely, 

V2.  _±d^ 

c2  Sr2  c,  dr 

where  //is  the  heat  flux.  A typical  transient  solution  of  this  equation  has  a form 

p(r.f)  = ■ X const,  x /(/), 

‘\  >2 

in  which  Eg  is  the  cosmic  particle  energy,  T\  is  essentially  the  time  of  heat  expansion,  and  Tj 
is  the  interval  between  the  peak  of  compression  and  the  "peak"  of  rarefaction.  When  the  gen- 
eralized model  applied  to  speciflc  cases,  the  differences  between  these  specific  models  is  found 
to  rest  on  choices  of  T\  and  72-  Actually,  EglT\  determines  the  magnitude  of  the  peak 
compression,  and  T2  determines  the  dominant  portion  of  the  frequency  spectrum  of  the  radiat- 
ed wave.  Measurements  in  laboratory  experiments  show  Tj  to  be  10  to  1 ns.  leading  to  the 
conclusion  that  the  significant  spectrum  of  p(r,i)  begins  somewhere  near  25  kHz  and  extends 
thereafter  upward. 

A crucial  question  is  the  range  limit  of  detectability  of  these  bursts  of  sound.  Local  noise, 
frequency  content  of  signal,  incident  cosmic -particle  energy,  propagation  attenuation,  etc.  all 
must  be  considered  in  determining  the  limit  of  "audibility".  Rough  approximates  have  been 
made  in  the  body  of  this  report.  For  example,  if  the  pressure  pulse  is  processed  in  a filter  cen- 
tered at  25  kHz  and  10  kHz  wide,  then  the  limit  of  range  in  units  of  meters  per  electron-volt  at 
which  the  transient  signal  magnitude  is  equal  to  the  local  noise  (also  at  25  kHz)  is  given  by  a 
particular  model  as  2 x 10  Thus,  in  this  model,  a total  cosmic  particle  energy  of  some 
5x10'“'  eV  is  needed  to  permit  detection  at  a limit  distance  of  100  m.  Other  models  yield 
somewhat  different  (or  vastly  different)  values,  depending  on  their  assumptions.  The  conclu- 
sion is  that  several  major  problems  in  modeling  remain  to  be  solved: 

• The  basic  physics  of  the  conversion  of  cosmic  particle  energy  into  acoustic  energy  en- 
visages the  formation  of  a shock  wave  due  to  nearly  instantaneous  rise  of  heat  deposition,  yet 
modeling  to  date  is  based  on  linear  acousatic  theory  of  small  amplitude  wave  motion. 

• In  all  mathematical  models  the  key  parameter  of  time  estimation  appears  explicitly, 
but  the  magnitudes  of  T\  and  T2  to  be  assigned  to  the  models  are  largely  guessed  at. 

• Experiments  conducted  to  date  verify  certain  parameter  dependencies  which  tend  to 
verify  particular  models.  These  dependencies  show  that  the  acoustic  pressure  is  proportional  to 
the  energy  of  the  particles  coming  in  from  space,  being  proportional  to  the  coefficient  of  ther- 
mal expansion  of  the  medium,  inversely  proportional  to  the  specific  heat  of  the  medium,  and 
inversely  proportional  to  the  distance  of  the  observation  point  (spherical  spreading).  In  addi- 
tion it  is  verified  that  the  radiation  is  sharply  pancake  shaped  normal  to  the  axis  of  the  cascade. 
Still  unknown  however  is  the  basic  mechanism  of  energy  conversion,  the  modeling  to  date  be- 
ing purely  phenomenological. 
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Fig.  2 - Addition  of  random  noise 
from  two  showers 


One  must  first  construct  an  acoustic  mode)  of  one  sfiower.  This  is  covered  in  the  main 
body  of  the  report.  Secondly  one  must  find  a way  to  add  sources  together.  In  view  of  the  ran- 
dom nature  of  the  showers  these  effects  must  be  added  on  a probabilistic  basis.  A schematic 
way  of  doing  this  is  as  follows.  Figure  2 shows  two  contributing  showers  radiating  nr  ise  to 
point  r.  Shower  1 is  the  jth  random  shower  of  the  kth  random  high-energy  particle,  and 
shower  2 is  the  1th  random  shower  of  the  mth  random  high-energy  particle.  Many  more  such 
showers  contribute.  The  total  time  record  of  reception  at  r from  all  contributing  showers  is  a 
simple  sum  of  time-delayed  elementary  pulses: 


F(r,/) 


r (i). 


Since  F is  a random  function  of  space  and  time,  we  will  attempt  to  construct  its  autocorrelation. 
First,  let  exactly  K contributing  particles  arrive  in  the  time  interval  O.T,  and  let  exactly  7 con- 
tributing showers  of  these  particles  occur  in  spatial  volume  V.  Then  the  total  pressure  due  to 
these  showers  is 


= Z Z Pjk 

j -1  k -1 


t — 


Ir  -R/a  ('a)I 


Now  let  g(J),  and  g(K)  be  the  probabilities  that  exactly  7 contributing  showers  will  occur  in 
volume  V and  exactly  K contributing  particles  will  occur  in  the  time  interval  O.T.  The  auto- 
correlation ill  of  the  received  time-pressure  record  at  point  r is  then 

■mi.t)  = i £ gu)g(K)  <fpjKiT.t)  pj^(t+\.i+t)  y . 

J-O  K-O  ' 

in  which  the  overbar  indicates  the  time  average  and  the  brackets  indicate  the  ensemble  aver- 
age. 
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in  which  ^(A")  is  the  probability  that  exactly  X events  will  occur  in  the  interval  0.(A  and  \ is 
the  average  number  of  events  per  unit  interval.  It  is  further  plausible  to  assume  that  space  and 
time  probabilities  are  statistically  independent.  Following  Rice  [2],  we  eonclude  that 

^>/i(l,7)  ^ f J p(i.i)p(T  +1,1  +r  ) dfdi  + ^ 

— oe  —op 

where 

. oo 

Pit,t)  / ^ = f f pix.t)  drdi. 

—oo 

The  space  and  time  averaged  power  spectrum  WiK.w)  is  then  obtained  by  Fourier 
transformation  of  ih: 

oo  oo 

^l+(K,<u)^  =s  J*  d-  J d I cos  (uT  cos  K I 

— oo  — oo 

= 4 yV^A'JSlK.to)!  2 + 4<(p7r7T)>  ^ 8(K,ru), 

in  which 

S{K,ai)  = J*  J"  piT,i)e'^^  ' drdi. 

— oo 

This  means  that  one  can  obtain  the  /jotver  spectrum  of  the  acoustic  noise  by  first  obtain- 
ing the  space-time  spectrum  5 of  the  transient  pressure  p and  the  space-time  average  of  the  to- 
tal pressure  and  then  adding  them  with  the  coefficients  as  shown.  Other  signal-processing 
schemes  can  be  used  to  obtain  further  information  on  total  effects  in  the  ocean. 


Models  of  the  Separate  Sources 

The  subject  of  this  report  is  the  development  of  acoustic  models  of  the  individual 
showers  generated  by  a high-energy  particle.  In  the  process  of  this  development  we  will  bring 
out  the  important  property  of  directionality  of  the  acoustic  signal  which  is  a prime  goal  of 
analysis. 

High  energy  atomic  particles  are  decelerated  upon  passing  through  sufficient  thicknesses  j 

of  solid  or  liquid  materials,  losing  kinetic  energy  as  they  do  work  on  adjacent  molecules  of  the  i 

medium.  The  energy  loss  is  assumed  to  take  place  in  the  following  steps:  A fast  moving  parti-  f 

cle  strips  off  an  electron  at  high  energy,  and  then  this  electron  strips  off  further  electrons  in  ' 

the  atoms  along  its  path,  the  process  continuing  in  a buildup  of  a multitude  of  fast-moving  ion- 
ized atoms  and  electrons.  These  collision  products  distribute  their  energies  among  other 
molecules  of  the  medium  by  additional  successive  collisions.  Each  loss  of  energy  reappears  ul- 
timitely  as  heat  in  the  absorbing  medium  or  as  a :hange  in  the  medium's  momentum  or  local 
chemical  state.  The  energy  transfer  occurs  quickly  and  is  thought  to  be  the  origin  of  experi- 
mentally observed  acoustic  transients. 

There  are  several  origins  of  high-energy  particles.  Man-made  sources  are  nuclear  reac- 
tors, particle  accelerators,  atomic  collisions,  etc.  In  the  natural  world  it  is  known  that  cosmic- 
particles  pervade  the  entire  universe.  Those  that  find  their  way  into  the  earth's  atmosphere, 
oceans,  and  crust  are  thought  to  originate  from  distant  galaxies,  supernovas,  and  the  sun. 

Their  detection  is  of  scientific  and  military  interest.  Several  schemes  of  detection  have  been 
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proposed,  depending  on  the  nature  of  the  particles.  If  the  particles  are  charged,  one  can  detect 
them  in  cloud  or  bubble  chambers  in  the  presence  of  magnetic  fields,  or  by  Cerenkov  radia- 
tion. If  the  particles  are  not  charged  (say  neutrinos),  their  interaction  with  other  matter  may 
induce  formation  of  secondary  particles  which  then  can  be  delected.  Finally,  it  has  been  pro- 
posed that  high-energy  particles  generate  acoustic  pulses  in  media  and  that  these  pulses  are 
sufficiently  intense  to  be  delected,  the  models  of  this  acoustic  pulse  generation  being  the  sub- 
ject of  this  report. 

Several  theoretical  models  have  been  proposed  to  explain  the  mechanisms  of  acoustic- 
pulse  formation  by  cosmic  particles.  The  dominant  models  are  thermoelastic  models;  in  these 
models  kinetic  energy  is  converted  into  heal,  which  raises  the  temperature  of  the  medium 
thereby  altering  the  medium  density  on  a short  time  scale  which  causes  a pressure  shock  wave 
to  radiate  outward.  Three  distinct  heat  models  have  achieved  prominence  among  scientists 
who  have  examined  the  problem  HI.  These  are  the  heated-rod  model,  the  heated-spot  model, 
and  the  heated-filament  model. 

In  the  heated-rod  model  the  cosmic  particle  is  thought  to  generate  a cascade  of  ionization 
products  that  instantly  heat  a cylindrical  rod  in  the  medium  of  length  L and  diameter  d (the 
numerical  values  of  L and  d depending  on  the  energy  of  the  particle)  to  a fractional  degree 
above  ambient  (the  temperature  rise  also  depending  on  particle  energy).  Typical  numbers  are 
L = 5 to  10  m,  d = 1 to  12  cm,  and  temperature  rise  = 1 for  the  case  of  a 10  eV  neutri- 
no. The  heated  cylinder  expands  in  10  to  0.01  fxs,  radiating  outward  a transient  shock  wave 
with  an  efficiency  of  energy  conversion  of  about  10  to  10  The  radiation  along  the  en- 
tire length  L is  coherent  over  the  range  of  the  bandwidth  frequencies  of  the  transient  and  is 
directional  with  angle  of  detection. 

In  the  heated-spot  model  the  kinetic  energy  of  the  particle  is  converted  explosively  into  a 
concentrated  heal  source  (the  hot  spot),  with  the  heat  then  diffusing  outward  in  finite  (but  ra- 
pid) lime.  This  thermal  shock  results  in  an  elastic  stress  wave  which  radiates  outward  as  a 
shock  wave  with  a speed  initially  greater  than,  but  eventually  equalling,  the  speed  of  sound  in 
the  medium.  The  pressure  pulse  of  these  waves  exhibit  compression  only  and  the  spectral  fre- 
quencies of  the  transient  propagate  independent  of  angle. 

In  the  heated-filament  model  the  healing  is  confined  to  narrow  filaments  (say  1 /ant  in  di- 
ameter), there  being  showers  of  them.  The  temperature  rise  in  a filament  (like  that  of  the 
hot-spot  model)  is  100  to  10,000  times  higher  than  that  in  a rod.  Each  filament  radiates  a 
separate  acoustic  transient,  and  it  is  assumed  that  all  the  filaments  in  a shower  radiate 
coherently.  The  heated-filament  model  (as  well  as  the  heated-spot  model)  predicts  acoustic 
transients  several  orders  of  magnitude  larger  than  the  heated-rod  model  because  of  the 
confinement  of  heating  to  small  volumes  (with  consequent  greater  temperature  rise)  and  be- 
cause of  the  much  shorter  time  of  heat  deposition. 

In  addition  to  these  three  thernioelaslic  models,  another  model,  not  readily  calculable  but 
of  equal  importance,  is  the  microbubble  model.  Here  the  concentrated  heal  of  a secondary 
product  particle  quickly  boils  the  liquid  into  a vapor  bubble.  The  growth  and  collapse  of  this 
bubble  radiates  sound  in  the  form  of  shock  waves. 

Thus  the  complexity  of  cosmic  particle  showers  makes  acoustic  modeling  highly  uncer- 
tain. The  most  satisfactory  approach  is  to  simplify  the  final  acoustic  radiator  to  a collection  of 
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recognizable  (and  computable)  acoustic  shapes.  In  this  survey  we  will  develop  in  detail 
cylinder,  sphere,  and  vapor-bubble  models,  with  the  cylinder  and  sphere  models  being  ther- 
moelastic models.  .These  models  will  be  used  to  make  numerical  calculations. 

Cosmic-Particle  Showers  in  the  Atmosphere  131 

The  dominant  high-energy  particle  arriving  from  outer  space  is  the  proton.  Approximate- 
ly 10 arrive  per  sec  on  a 1-km^  area  at  the  top  of  the  atmosphere  with  energies  equal  to  or 
greater  than  10^  eV.  Only  the  higher  energy  protons  will  reach  sea  level.  The  energy  flux 
spectrum  of  these  high  energy  protons  on  the  top  at  the  top  of  the  atmosphere  is  extensively 
measured.  A few  entries  are: 


Energy  (eV) 

Flux  (s 

10'* 

10 

10'^ 

10 

10*9 

10 

Upon  entering  the  atmosphere  the  proton  collides  with  air  molecules  about  once  every  7 to  8% 
of  atmosphere  thickness.  Protons  with  energy  of  = 10*^  eV  create  about  10  pi  mesons  at  col- 
lision, and  those  with  energy  of  = 10 eV  about  20  pi  mesons.  A proton  survives  each  colli- 
sion with  about  half  its  initial  energy.  Pi  mesons  themselves  undergo  collisions,  creating  addi- 
tional pi  mesons  in  the  same  numbers  as  proton  collisions.  Thus  a cascade  of  pi  mesons  is 
created  by  each  initial  proton  collision,  the  number  sequence  being  one  proton,  then  ~ 10  pi 
mesons,  ~ 100  pi  mesons,  ■—  1000  pi  mesons,  etc.  Since  the  atmosphere  is  relatively  thick, 
most  of  the  proton  energy  that  goes  into  creating  pi  mesons  is  absorbed  before  sea  level  is 
reached.  However  the  newly  created  pi  mesons  are  unstable  and  immediately  decay.  The 
plus-  and  minus-charged  pi  mesons  decay  into  plus-  and  minus-charged  muons  and  neutrinos. 
The  neutral  pi  meson  decays  into  two  gamma  rays.  Decay  products  then  interact  with  the  at- 
mosphere to  produce  a cascade  of  electrons  and  protons.  Muons  do  not  interact  strongly  as  do 
pi  mesons  and  protons,  so  they  penetrate  to  sea  level  and  below.  Neutrinos  have  even  weaker 
probability  of  interaction  with  the  atmosphere.  They  too  penetrate  into  the  ocean. 

If  a high  energy  neutrino  (say  10 '“’eV)  does  interact  with  the  atmosphere,  about  half  the 
initial  energy  goes  into  a hadronic  cascade  (a  cascade  of  "hard”  particles  possessing  significant 
mass),  the  other  half  going  into  the  creation  of  muons  or  electrons.  The  hadronic  cascade  has 
the  following  history:  The  neutrino  creates  pi  mesons.  Of  the  total  population  of  charged  me- 
sons within  50  cm  along  the  interaction  path  about  10  generate  (say)  10  more,  resulting  in  a 
shower  of  mesons.  Pi  mesons  with  zero  charge  create  a positron-electron  cascade  spatially  dis- 
tributed in  a cone  estimated  to  be  5 to  10  cm  across  and  5-10  m long.  The  number  of  particles 
in  this  cone  exceeds  10'*. 

If  a high-energy  muon  (say  10 eV)  interacts  with  the  atmosphere,  it  creates  a particle 
shower  approximately  every  100  m.  These  showers  may  begin  as  gamma  rays,  positron- 
electron  pairs,  or  pi  mesons.  The  energy  in  a shower  is  about  10  ” eV,  there  being  about  100 
particles  in  a cone  a few  meters  long  and  about  10  cm  wide.  If  the  high-energy  muon  has  10 
times  the  preceding  energy  (~  10'“'  eV),  it  will  also  generate  showers  every  100  m but  which 
will  be  10  times  more  energetic  with  10  times  more  particles.  Along  each  100-m  segment  of 
the  shower  path  there  is  a regular  ionization  loss  of  2 x 10*  eV  per  cm,  or  2 x lO'^eV.  This 
is  roughly  1%  of  the  initial  muon  energy  and  can  be  assigned  as  a 1%  energy  feed  into  each 
shower. 
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Statistics  of  Particle  Showers  Reaching  Sea  Level 

A useful  model  of  an  air  shower  reaching  sea  level  is  a thick  rod  of  area  10  m by  10  m, 
one  rod  occurring  per  square  kilometer.  The  flux  of  particles  in  this  rod  is  a random  quantity 
but  can  be  modeled  as  follows: 

• In  every  100  s there  will  occur  a shower  of  10^  particles.  Of  these  some  99%  are  posi- 
trons or  electrons,  which  are  stopped  in  approximately  10  m of  water.  About  1%  are  muons, 
which  penetrate  deeper. 

• In  every  3 hours  there  will  occur  a shower  of  10*  particles. 

• In  every  couple  of  months  there  will  occur  a shower  of  10 particles.  Of  these  some 
10*  are  muons,  10 of  which  have  enough  energy  to  penetrate  5 km  into  the  ocean. 

Each  particle  has  a mean  energy  of  10^  eV.  The  particles  which  penetrate  into  the  ocean 
deep  enough  to  be  of  significance  in  this  survey  are  the  neutrinos  and  the  muons.  The  flux- 
and-energy  statistics  for  the  neutrinos  are  briefly  as  follows: 

• One  per  year  interacts  in  1 km  ^ of  water  with  an  energy  of  10  '•  eV, 

• 100  per  year  interact  in  1 km  ^ of  10  eV,  and 

• 10“*  per  year  water  with  an  energy  of  10 eV. 

The  flux-and-energy  statistics  for  the  muons  are  briefly  as  follows: 

• 300  per  second  arrive  on  1 km  ^ of  ocean  at  10 eV,  each  accompanied  by  100  muons 
at  10’ eV, 

• One  per  second  arrives  on  1 km^  of  ocean  at  10 eV,  accompanied  by  ~ 1000  mesons 
at  10’ to  10'‘’eV. 

• One  per  100  s arrives  on  1 km  ^ of  water  at  10  eV. 


THERMOELASTIC  MODELS 

Sound  Sources  and  Their  Governing  Equations 

The  generation  of  sound  by  cosmic  particles  in  the  ocean  may  be  investigated  theoretical- 
ly by  adopting  a reasonable  model.  Since  the  mechanism  of  energy  deposition  and  acoustic 
conversion  is  not  known  with  certainty,  it  will  be  useful  to  review  various  potential  models. 

A first  approach  to  modeling  is  to  state  the  equation  of  state  for  liquid  pressure  /?in  which 
small  deviations  from  the  ambient  (or  equilibrium)  pressure  />„  are  dependent  on  changes  in 
density  Ap  and  in  extropy  AS! 
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For  a heat-conducting  viscous  fluid  the  excess  pressure  (acoustic  pressure)  is 

' 2 I * 
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in  which  c is  the  speed  of  sound,  tj  and  tj  are  viscosity  coefficients,  x is  the  thermal  conduc- 
tivity (units:  N/K  ■ s),  and  Cy  are  specific  heats  at  constant  pressure  and  constants  volume 
(units:  ■ 5^),  and  v is  the  vector  fluid  velocity.  The  simplest  model  of  acoustic  pressure 

generation  is  the  one-term  statement 


P 


= C'V  • 


This  is  the  case  of  a linear  nonviscous  fluid  in  which  the  conduction  of  heat  is  negligible  and 
the  mechanism  of  sound  generation  is  exclusively  a function  of  change  in  density.  However, 
the  generation  of  sound  pulses  by  cosmic  particles  is  a complicated  phenomenon.  Hence  a 
more  complicated  model  will  be  needed. 

We  require  various  more  complex  mechanisms  for  changing  density  in  liquids.  These 
mechanisms  are  customarily  listed  as  "sources"  in  writing  the  general  linear  equation  of  the 
propagation  of  sound  (4,  p.324l: 


V^p  — f)K  = — /'(r,  t)  -h  div  F — V 


T ■ V + 


in  which  thermal  and  viscous  terms  have  been  omitted  on  the  left-hand  side  but  are  included 
on  the  right-hand  side  as  potential  sources.  Here  /(r,  1)  is  the  sum  of  all  monopole  sources, 
div  F represents  the  dipole  source  terms,  and  V ■ F • V is  the  quadrupole  source.  Our  interest 
will  be  focused  on  the  monopole  term.  The  general  formula  for  the  monopole  sources  is  given 
by  Eq.  7.1.22  of  Ref.  4; 


aypK  be 

bx 

bP  oye 

br  '1 

df 

br 

Sk 

bP 

aye 

2 

bP 

T 

bt 

di  ' 

_iL 

y 

Here  q is  the  volume  flow  (units:  m^/s  ■ m^)  and  P \s  the  equilibrium  pressure.  Other  symbols 
will  be  explained  below.  The  first  monopole  term  represents  an  introduction  of  new  fluid  into 
the  medium,  the  second  term  represents  an  introduction  of  time  varying  heat,  and  the  third 
and  fourth  terms  are  sources  by  virtue  of  the  temporal  and  spatial  nonuniformity  the  medium 
properties.  Although  other  sources  are  possible,  most  modelists  of  the  acoustic  effects  of  cos- 
mic particles  in  liquids  consider  the  heat  source  as  the  most  appropriate: 


/(r,  i)  = ayK  ~ 
dr 


III- 
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in  which 


dT 


(units: 


N 


nr 


yx 


(units:  m^/N);  y =—f-. 


C„ 


dP 


C. 


Since  € is  the  heat  added  per  unit  mass  per  unit  time,  the  factor  f)t  is  the  heat  added  per  unit 
volume  per  unit  time.  The  factor  e/C^  is  the  temperature  change  due  to  the  addition  of  heat. 
Introducing  the  coefficient  of  thermal  expansion  fi  ( units:  K one  can  write 

= li- 


ayK  = 


llL 

dr 


Thus,  by  regrouping  all  factors,  one  arrives  at  a more  significant  expression  for  the  monopole 
source  directly  attributable  to  heat: 


...  , H BH 

H = fj€ 


9f 


units: 


N ■ m 


;ir 


The  physical  significance  of  the  source  is  this:  The  addition  of  heat  results  in  a change  in  tem- 
perature, which  in  turn  causes  an  expansion  of  fluid,  thus  altering  density;  and  the  change  in 
density  generates  an  acoustic  pressure,  as  indicated  by  the  first  equation  of  this  section.  The 
equation  of  propagation  of  the  acoustic  pressure  is  therefore 


di^ 


dt 


The  negative  sign  appears  on  the  right-hand  side  to  insure  a positive-sign  source  when  the 
temperature  increases  {ft  being  negative). 


It  will  be  useful  to  picture  the  acoustic  phenomenon  described  by  this  equation.  In  anal- 
ogy with  a siren  which  ejects  puffs  of  air  by  chopping  a steady  air  stream,  the  equation 
describes  a "heat  siren"  which  ejects  puffs  of  heat  by  chopping  a steady  deposition  of  heat.  The 
heat  stream  is  H (calories/s  per  unit  volume),  and  the  chopping  action  is  d/dr.  If  the  chopping 
is  steady  state,  one  sets  d/dt  — —io>.  A "one-shot"  siren  generates  a transient  heat  puff. 
These  heat  puffs  in  turn  generate  density-temperature  puffs  (through  multiplication  by  C^”'), 
and  the  latter  are  converted  to  pressure  (density)  puff's  by  ft.  The  intensity  of  the  sound  gen- 
erated will  be  directly  proportional  to  the  magnitude  of  the  steady  heat  flu.x  (=  H)  and  to  the 
rapidity  of  chopping  ( = d/dr). 


This  picture  is  common  to  all  ihermoelastic  models.  As  we  shall  see,  a different  model  is 
constructed  for  each  description  of  H and  d/dr.  Actually  H is  specified  as  the  calories  per 
second  per  unir  volume,  making  the  developed  pressure  obtainable  only  by  integrating  over  the 
effective  volume  which  carries  the  heat.  The  heat  in  question  is  of  course  not  directly  deposit- 
ed by  the  cosmic  particle,  which  carries  kinetic  energy  only.  However  for  making  estimates 
the  magnitudes  of  kinetic  energy  and  the  developed  heat  are  made  equal. 

We  will  now  survey  several  useful  models  and  will  provide  development  of  details  when 
it  is  judged  useful  to  an  understanding  of  the  limits  of  validity. 
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Heat  Models  Based  on  — (l/c^)  (6 ^ = —(jilCp)  (dH/di) 

We  first  enquire  into  the  nature  of  H.  Let  the  initial  energy  of  a cosmic  particle  be 
(units:  N/m).  We  assume  that  when  the  particle  enters  the  ocean,  the  heating  ^effect  on  the 
water  becomes  volumetric.  We  can  represent  this  effect  as  where  £„  is  the  state- 

ment of  the  energy  disturbance  upon  the  water,  as  distinct  from  the  particle  energy  Next 
we  assume  that  at  any  particular  point  r„  along  the  track  the  water  is  at  first  quiescent,  then 
responds  to  the  heat  disturbance,  and  then  is  quiescent  again.  The  heat  description  is 


A//  = pulse  X 


A£o 


(units:  ,\m/s  ■ Three  specific  forms  of  this  function  will  illuminate  the  meaning  of  the 
symbols. 

For  a point  (or  delta)  deposition  in  space  and  time,  the  form  is 

A£' 


A//  = pulse  X 


AT,. 


= £„  8(r„)  S(/  ). 


For  a deposition  that  is  delta  in  time  and  uniform  in  cross-sectional  area  S,  the  form  is 

1 clE(=„) 


S dz„ 


8 (L. ) 


1 


which,  if  £„  (r„)  = £„  e becomes 


AW  = 


e "'yU,,)  (easel) 


and,  if  £„  (r„ ) is  Gaussian  over  volume  becomes 

{.-  -£72)^8 


AW  = — expl 


8(t„).  (casein 


For  a deposition  that  is  delta  in  time  and  varies  with  and  with  S(ri),  rj.  = V.v^  + , the 

form  is 

dE„ 

“ ' H 

which,  if  the  variation  with  ris  parabolic,  becomes 


AW  = 


1 - 


8(f„)  (casein) 


and,  if  the  variation  with  ris  Gaussian,  becomes 


dFAz„) 

AW  = — ; — exp 


dA, 


(case  IV ) 
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In  these  three  forms  the  pulse  nature  of  A// is  written  as  S(f„).  Other  pulse  expressions  may 
be  more  appropriate.  Writing  "pulse"  as  hii)  ( units:s  we  list  a few  pertinent  forms: 


h,  (t)  = rect  (0,  T)  or  rect  (-T/2,  T/2), 


“triangle  (0,  T)  or  triangle  (-  T/2,  T/2), 

/ijit)  =cos  11 /rect  (- T/2,  T/2). 

Each  pulse  function  has  a Fourier  spectrum  of  frequencies.  Letting  u and  s be  Fourier  pairs, 
we  have 


rect  u — sine  s. 


triangle  u -*  sine  ^ s, 
cos  TT 


rect  u — > -r  sine 

'+1 

-t-  — sme 

s 

2 

2 

2 

2 

Since  time  is  the  analog  of  u and  frequency  is  the  analog  of  s,  pulse  descriptions  and  frequency 
descriptions  are  interchangeable  by  means  of  Fourier  transformation.  It  will  be  convenient  in 
the  following  analysis  to  consider  both  temporal  and  frequency  descriptions.  In  all  cases  where 
the  steady  state  solution  p(i|<u)  of  the  pressure  has  been  calculated,  the  temporal  solution  will 
be  the  standard  one. 


p(x,  t)  = f p(t\oi)  h{o))  e 

''  -<«  2tt 


in  which  h(<u)  is  the  Fourier  transform  of  the  pulse  h(t). 
We  will  now  survey  useful  models. 


Let  the  initial  flux  of  particles  entering  the  ocean  be  a beam  of  intensity  /„  (units: 


Ntn 


-I 


m ^).  The  word  initial  means  that  we  know  /„  at  some  particular  field  point  be- 


fore additional  absorption  takes  place.  We  assume  that  the  rate  of  absorption  with  distance  is 
exponential  of  form  e where  « (units:  m “')  is  the  absorption  coefficient.  This  space  rale 
of  absorption  is  precisely  H,  namely, 

-d 
dz 

The  time  variation  of  H is  as  yet  unspecified.  Initially  we  resolve  the  time  variation  into  har- 
monic components.  Thus  the  complete  specification  of  the  energy  deposition  in  the  fluid  be- 
comes 


H = 


■"-) 


W = « /„  c 


Next  we  must  assume  a geometrical  description  of  the  space  in  which  the  acoustic  field  will  be 
generated.  For  our  first  model,  which  we  call  heat  model  1,  the  space  is  assumed  to  be  infinite 
in  all  directions.  We  will  develop  this  model  in  detail. 


Heat  Model  la  On  Unbounded  Space)  iWestervelt-Larsen-Hanishi 


The  transient  problem  to  be  solved  is 

i 1 /) 

V^p(r,i) r — - — ; — = 


-<7(r„. /fl).  q 


HH  fi 


Hr  C.. 


1 


13 


S.  HANISH 


The  solution  of 


is  known  to  be 


V"/> 


2„  _ _L 


— — -r  = -4nq  . q = ql^-n . 
tr  dr 


/ 

4ttp(t.  I)  ^ di„  J dV„  G'(r,  r„lf.  t„)q'(r„.  i). 

where  G is  the  Green’s  function  of  the  scalar  wave  equation  (5.  p.  834].  For  unbounded  space 

h[R/c  - U - / .)] 


G = 


where 


and 


= (.V  - .v„ ) ^ + Cv 

Thus  the  solution  of  the  transient  problem  is 

dV 


R 

R =\t  - r„| 

>•» ) ^ + (--  - 


r r li  C'^  . r,,SlR/C-U-,„)]dH{T„.rJ 


477- C„ 


J dy„ 


R 

- R/c) 

R 


dr 


in  which  the  prime  signifies  differentiation  with  respect  to  the  argument  of  H (namely 
t — RIO.  Let  the  volume  of  integration  be  a cylinder  of  area  S,  length  /.  Then,  substituting 
the  previously  derived  form  of  H in  which  the  transient  state,  is  replaced  by  the  steady  state 
we  arrive  at 

ic^  a l„  Sp  dz„  e e * 

= 7—7^ J 7 • 


47r  r„ 


In  the  far  field 


/?  = I r|  — cos 

where  |r|  is  the  distance  from  the  origin  to  the  field  point.  After  a certain  length  /*  = n 

is  effectively  negligible  in  its  contributions  to  the  value  of  the  in- 


the  function  e 


^ -1 


tegral.  By  making  / > / , one  can  extend  / to  infinity  and  thus  approximate  the  integral  to  read 

,-A|r|  P 


p(r,<o)  = 


—oja  Ig  S P . , , « -loii 


477  Cp  |r| 


d + ik  cos  P ' 

The  magnitude  of  pressure  has  a directivity  F(«),  where 


k =^. 


, , fo  ■5/3 

I/"!  = 7 -1  r Rio).  F(P)  = 


1 


477  C^|r| 


vr 


+ l*^cos  ^ a 


The  maximum  occurs  at  broadside  (P  =90°).  (See  Ref.  6 additional  details.)  Although  the  ex- 
ponential heat  depostion  is  plausible,  other  types  of  track-dependent  deposition  are  equally  so. 
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An  example  is  the  Gaussian  dependence  on  penetration  length  (=  c„),  which  earlier  was  la- 
beled case  11.  Thus  we  use  the  earlier  equation 

- Ul ) ^8 


AW  = — exp 


Z.2 


and  consider  only  a single  frequency.  Defining  the  intensity  flux  as  l„(r^)  (energy  flux  per 
unit  cross  sectional  area),  we  write 


\H  * — exp 


-Z./2)28 


„ -lOJi 


This  implies  that  the  "Gaussian  length"  of  deposition  is  L,  that  deposition  is  a maximum  at  Z./2, 
and  that  the  standard  deviation  of  deposition  is  tr  = ^(1/2)  (Z.^/8)  = Z./4.  The  radiated 
acoustic  pressure  is 

I 2g 


/’(r.  <«)  = 2 "Tr'  f 

4rr£  C„ 


exp 


- Ul ) ^8 


exp 


-H  I ) 


l„  (r,  ) <n  \, 


For  simplicity  we  take  to  be  a constant  independent  of  the  cross  sectional  area  variable 

r^.  Thus,  dy  = dz„  dS(r„  ),  and 


, , n lo  s 

p(t,o>)  =~ — — -e 


477  L C„ 


L exp 
-wj 


(r„  - Ul ) ^8 


Ul 


,ikH 


cfz  k — (jj/c. 


0 


In  the  far  field 


/?  ==  |r|  - r cos  0, 


where  « is  the  angle  of  the  observation  vector  with  the  ; (or  track)  axis.  Let  ,v  = z^,  - Ul,  dy  = 
dz^,-,  then,  extending  the  limit  of  integration  lo  infinity. 


-i(u  fi  /„  S 


P(t,  w)  = ^ ^ e -dk/DI-cosH  J ^ -y^/U  -ikycmH 


477  Z-  |rl 


Let  8/Z.^  = a and  ik  cos  P = f3.  Then  the  integral  reduces 


to 


oo 

C ^ —ti  (.F  ) 

'\j 

0 

2a 

= ,,(» (/i/2..  v'tf 

2v;7 


_ L [7 

4 V 2 ‘ 


.12 


From  this  the  magnitude  of  pressure  is 

lp(r.to)l  ~ ^ ^ 


lr| 


2770 


In  most  cases  LA  >>  I,  so  that  the  directionality  of  the  pressure  is  sharp  in  the  broadside 
direction  («  ==  90°).  The  quantity  /„Sis  again  equal  to  E,JT, , the  "time"  7'.  being  an  estimate. 
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i 

I 


I 

j 

[ 


I 


1 


A second  example  of  a spatial  deposition  of  energy  which  is  other  than  exponential  is  the 
uniform  deposition  in  which  dEJdz  is  constant.  Thus 


so  that 


In  the  far  field 


or 


“ T — «(/<,) — — e 

•S  dZg  S dZg 


pit,  tu) 


-m  <3  1 dE  fe^ 

4n  Cp  r,  dz„ 


I R 


IH  TtI 

Q h r J 


ai  f)  (dEg/dz) 

1^1  ^ T I I ^ 

477  Cp  r,  |r| 


n L 

sin 

cos  0 

X 

77  Z. 


COS  9 


Here  again  is  unknown.  The  radiation  pattern  is  sharp  at  broadside,  just  as  in  the  case  of 
the  Gaussian  deposition. 


Heat  Model  Ib  Gn  Unbounded  Space) 

^ When  the  heat  deposition  is  transient  and  there  are  no  preferred  frequencies,  we  replace 
e by  "pulse"  h(i).  The  solution  is 

„(r  ,)  =_J_  f ^ 

477  Cp  R 

Since  all  physically  compatible  transients  hit)  start  and  end  in  finite  time,  the  time  derivative 
will  have  both  positive  and  negative  values  (/?'  = dh/dU  - Rfc)).  Hence  the  sound  radiated 
will  exhibit  both  compression  and  rarifaction. 

A simple  application  of  this  formula  is  to  assume  the  heat  deposition  is  a delta  function 
in  space  and  time; 


~ ^0  ^ S (•'o ) units: 


N • m 


s ■ m 


3 


Then  the  radiated  transient  is 


P^T,  t)  = 


4’^  Cp  -0 


5 d'o  f dy„ 


I A , V 

7 - 


IrL 

477  C„ 


8' 


R 

1^-, 


8'(/,)8(r.) 
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Eig.  3a  - A double  impulse  of  pressure 


This  pressure  history  exhibits  a double  impulse  at  point  |r|  (Fig.  3a).  The  first  or  compressure 
pulse  is  physically  due  to  the  shock  of  expanding  heated  fluid  acting  against  the  inertia  of  the 
stationary  non-heated  fluid;  the  second  or  rarifaction  pulse  is  due  to  the  "vacuum  wave  creat- 
ed by  the  sudden  cessation  of  expansion  of  the  heated  fluid  while  the  inertia  of  the  non-heated 
fluid  continues  its  outward  motion  (Ref.  7,8).  As  a second  example  we  assume  a linw  history 
of  deposition  to  be  a sinusoidal  pulse: 


h(i) 


= (sin  ll/)n 


I - T/2 
T 


Here  II  is  a rectangle  function  of  unit  height  and  of  base  T s centered  at  772  s,  which  is  zero 
outside  the  range  0 to  T.  The  radiated  pressure  for  exponential  deposition  of  heat  is 


p(T.l)  = 


H Iq 
477  C„ 


S “K  - 


If  cos  If  (f  - )!”] 


(t  - RIc)  - Til 


-t-  sin  If  (i  — R/c)n 


0 - R/c)  - 772 


The  second  term  on  the  right-hand  side  represents  a pair  of  impulses  generated  by  the  step 
functions  appearing  inf”!.  We  shall  neglect  them  here,  but  will  consider  them  later.  Setting 
the  limits  of  integration  to  infinity,  we  find  the  pressure  at  distance  R > > /"/X  to  be 


p(r,  /)  = 


fi  If  « /„  s 

j«  cos  If 

— K cos  sin  11 

7-1^ 

477  (a^  -F  cos^  ff)|r|  1 

c 

c 

xn 


(i  - R/c)  - 772 
T 


K = 


c 


At  w=  90°  the  shape  of  the  impulse  pressure  is  shown  in  Fig.  3b.  We  see  again  that  the  pres- 
sure is  a maximum  at  n =90°  (broadside).  This  model  is  analogous  to  a "one-shoi"  siren.  It 
will  be  used  in  a later  section  on  numerical  calculations. 


Heat  Model  II  Qn  Semi-Infinite  Space) 


In  heat  model  II  one  again  assumes  that  the  cylindrical  shower  in  the  liquid  due  to  the 
interaction  products  of  a cosmic  particle  is  caused  by  an  external  beam  of  the  same  cross- 
sectional  area.  This  time  however  the  beam  enters  a semi-infinite  space  (z  > 0)  from  a 


l! 
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I'ig.  4 - Laser  beam  entering  an 
absorbing  half  space 


nonabsorbing  medium  (region  I)  into  an  absorbing  medium  (region  11)  (Fig.  4).  The  energy 
deposition  rate  is  again  taken  as 

H = a l„e 


in  which  the  units  of  H are  again  N • m/m^  • s and  the  units  of  /„  are  N • m/m^  ■ s.  Each 
point  in  region  11  along  the  centerline  of  the  beam  is  an  elementary  point  source  of  spherical 
waves.  Corresponding  to  it  is  an  image  point  source  in  region  1.  The  sum  of  these  two  waves 
integrated  over  the  length  of  the  real  beam  gives  the  field  pressure  at  r.  Using  the  basic  equa- 
tions of  model  1,  one  finds  the  solution  to  be 


p(r) 


— 13  a Ig  S i w 
47r  C„ 


; 


r -r 


|r  -^2! 


±. 


For  an  acoustically  soft  interface  one  chooses  the  negative  sign.  In  the  far  field,  where  R > > 
/‘A  (/  being  the  length  of  the  beam,  taken  as  a “'),  one  can  write 

I r - r|  I = /?  — z cos  9 and  | r - r2 1 = /?  + z cos  9. 

To  simplify  the  analysis,  the  integral  over  z is  extended  to  infinity.  Thus 
~(u  (3  a Ig  S gikR  -Kill 


p = 


477  C„ 


X 


-A;cos«  _ „ik: COSH 


1 dz 


or 

^ ~ ^ Ig  S p ikH  -Kill  2 kl  cos  9 

^ 477  R \ \ cos^  »■ 

(Kozyaev  and  Naugol’  nykh  ref.  7).  The  amplitude  of  radiated  pressure  is  then 


IpI  = 


13  Ig  S 
477  Cp  R 


F{9), 


F(9) 


2 kl cos  9 
1 -I-  (A7)  ^ cos  ^ 9 


The  significance  of  this  result  is  now  examined.  First,  the  amplitude  depends  directly  on  the 
power  of  the  beam  IgS,  where  5 is  the  cross-sectional  area.  This  is  actually  the  energy  deposi- 
tion rate  in  units  of  N m/s.  Second,  the  amplitude  depends  on  the  frequency  <0,  which,  in  the 
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time  domain,  represents  b/bi.  This  is  the  rate  at  which  the  energy,  once  deposited,  is  released 
to  the  medium.  Third,  the  amplitude  varies  inversely  as  /?,  showing  that  the  far-field  radiation 
is  spherical.  Although  the  beam  itself  is  a cylinder  and  the  near  field  is  cylindrical,  the  finite 
length  of  the  radiation  eventually  appears  as  a point  source  to  a distant  observer.  Fourth,  the 
pattern  function  is  a strong  function  of  angle  W:  at  M =90°  , (in  the  plane  r =0)  /■'(90°)  =0,as 
required  for  a pressure-release  interface.  The  pressure  reaches  a maximum  at  an  angle  H such 
that 


cos  « =77.  or  = cos  '((»/A  ). 
kl 


If  kl  is  large  (//A  > > 1),  the  angle  H lies  near  90°: 


0 ^ ^ - -L. 


Thus  the  maximum  radiation  occurs  near  90°  (horizontal  direction).  At  90°  the  radiation  van- 
ishes. The  amplitude  of  pressure  in  the  direction  of  the  maximum  is  independent  of  kl.  In 
contrast,  along  the  beam  in  the  direction  of  the  z axis  the  pattern  amplitude  is 

FiO)  = — — . 

I -f-  ikl)^  A/  A 

showing  that  the  pressure  amplitude  for  //A  > > I is  determined  by  the  absorption  a. 


A parameter  of  importance  in  the  detection  of  cosmic-particle  noise  is  the  total  power  ra- 
diated (an  integrated  elfect): 

(F  = 2-rrR^  f sin  0 cl  H 

i 2pC 


1 ‘0^  I . , , A/ 

= 7 y.  y. 77  arctan  A/  - . 

SrrpC  A/  1 + (a/)2 

Again  if //a  > > 1 , this  reduces  to  the  simple  form 

W — ; — Oi  <»:. 

16  (>  C'l 

The  halfwidth  of  the  main  lobe  of  field  intensity  is  then 

A7  A 

Until  this  point  one  has  assumed  that  the  incident  beam's  effective  cross-sectional  diameter  is 
small  in  compari.son  with  the  acoustic  wavelength  generated.  If  the  two  are  comparable,  there 
is  an  additional  directivity  arising  from  diffraction.  Thus  in  place  of  Fill)  one  has 

Fin.d,)  =F(H  )<).(«), 


</)(«)  = 


2 J\  (ka  sin  H ) 
ka  sin  fi 


The  maximum  value  of  this  function  is  unity  (at  H = 0).  Hence  the  diffraction  due  to  large 
beam  diameters  reduces  the  amplitude  in  all  other  directions. 
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Fig.  5 - Pressure  patterns  of  a laser  beam  in  an  absorbing  half  space 

A sketch  of  F(9)  (no  diameter  — diffraction  effects)  is  shown  in  Fig.  5 for  k/ small  and  kl 
large.  The  directionality  closely  resembles  a dipole  in  which  the  beamwidths  depend  on  kl. 

It  is  instructive  to  compare  these  results  with  the  steady  state  model  1.  The  average  in- 
tensity of  sound  generated  over  one  period  277/01  of  model  I is 

. ^ 1 4,  = 1 ^ Ip  ^ ^ 

(277/to)  -S  2 i>  c ‘ Ific  A ttRCj,  + k"^  cos"^  9 

Assume  the  effective  depth  of  penetration  is  / = a then  the  directivity  factor  F(9)  be- 


F(9)  = ^ . 

1 -t-  k^!^  cos  9 

The  maximum  intensity  of  sound  occurs  at  9 =90°  (broadside).  The  angle  at  which  the  inten- 
sity is  reduced  to  1/2  of  its  maximum  is 

H = JL  + iL 
«l/2  - 2 + A - 


The  total  power  generated  is 

tr  to  /a  2 .9,2 

W = — [arctan  kl  — arctan  ( —kl)  ], 

16  77  o C2 

in  which  for  A/  > > 1 (//A  > > 1 ) the  bracket  has  value  tt,  so  that 

» wfi^  (/„  9)2 

If  = r . 

16  ,,  C-2 

A sketch  of  F(9)  for  the  two  cases  A7  small  and  kl  large  is  shown  in  Fig.  6. 


Heat  Model  III  (Heated-Rod  Model) 

In  heat  model  III  the  entire  shower  of  interaction  products  from  one  cosmic  particle  is 
considered  to  be  an  elastic  rod  L cm  long  and  a cm  in  radius.  A mathematical  model  of  such  a 
rod  may  be  constructed  as  follows:  The  equation  to  be  solved  is 

_2  • 61 2;)  li  1 
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kl  small 

I'ig.  6 Intensity  patterns  of  a laser  beam 
in  infinite  space 

in  which  Q is  the  heat  flux  (units:  N m/s  m ^),  with  the  symbol  Q replacing  //  of  the  previous 
sections  to  avoid  confusion  with  the  symbol  ’ for  Hankel  functions.  We  shall  be  interest- 
ed in  time  harmonic  solutions  /„)  = q(r,,)e  so  that  the  above  equation  reduces  to 
the  inhomogeneous  Helmholtz  equation: 

+ k^p  = -4iTq(r„),  q{T„)  = - 

The  solution  is  easily  obtained  by  use  of  the  appropriate  Green’s  function  G: 

=Jj  (rlr„)  r/f„. 

in  which  is  the  volume  of  the  source.  Because  we  will  employ  cylindrical  coordinates,  it  is 
convenient  to  replace  both  pand  Gby  their  Fourier  transforms,  pir^.n)  and  FlGl  , 
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Here,  W = 0 is  broadside  and  m //,,*'  ’ (k*P)  is  the  two-dimensional  Green's  function.  Thus 


We  next  assume  the  heat  deposition  Q is  uniform  in  r over  the  length  L of  the  rod  and  zero 
outside: 


qU„)  =C»<r„^)rect  “ y- y 


Now 


/ 


l.  2 


-L2 


-vz» 

e rect 


_L 

2 ■ 2 


T sin  ITTn)  , nL 

■ V(>,)  ■ . 


kL  sin  H 
2 


We  must  next  evaluate  the  expression, 

(k'P)Q(>-JdS(r,J 

Only  the  asymptotic  field  will  be  of  interest. 

lim— °o  //")  (k*p)  = /— 

k'P  V TT 

In  the  far  field 


rAV 


■’.A*  = x/a ^ — T)^  and  Z'  = | — /•„ 


P-  r, 


r,,^  cos  (b, 


where  <t>  is  an  azimuthal  angle.  The  absolute  value  of  the  pressure  field  is  then 


where 


^ ^ ^ ■\/ irk\  I 

We  next  recover  the  field  in  cby  taking  the  inverse  transform  in  r: 

LA(r^.k*)e'^>=-^. 

Ir  ^ 2tt 

Since  the  first  zeros  of  sin  IV/ (V  are  tj  = ±1ti/L,  we  simplify  by  taking  A to  be  constant 
between  these  limits  with  a value  of  A (r^  , A)  (i.e.  tj  = 0 );  thus 

p(r  0)  = f LAir^.k)-^^  =2A{r  k). 

T 

— r ! f>  M.-  • 1,1 

ox  ‘ 


A (r^  ,k)  =77 


Q^r,„)  dSir) 


Noting  that  is  constant  across  the  endcap  areas  of  the  cylinder,  we  approximate  the  in- 

tegral (and  the  absolute  value  of  it)  by  QS.  From  the  definition  of  Q vie  then  multiply  n by 
L/L  to  obtain 

0 , .iiL 

477C„  47rC,,  477C„  k 
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SO  that 


’’  277  -\/  C^L 


The  choice  of  A is  somewhat  arbitrary.  Let  us  restrict  attention  to  frequencies  25  kHz  anci  less, 
so  that  the  smallest  wavelength  K = 1 .5  x 10^/25  x 10^  — 0.06  m.  Since  the  length  L of  the 
shower  is  about  10  m,  it  is  seen  that  A = L/160.  This  choice  is  equivalent  to  stating  that  the 
radius  of  the  first  Fesnal  zone  at  a distance  equal  to  L is  taken  to  be  Rf  = (aL  ) ’ ^ = L/47r.  or 
A =L/1677  ^.  Thus  the  field  will  be  measured  in  the  Fesnel  diffraction  zone.  Hence 

1 / L w2E\(o\li  1 f2Ef\ti 


p(r^.o)  = 


\6n^r, 


Cp^JLr\ 


Finally, 


Pir^.T))  = 


sin  fF 


In  the  DUMAND  Proceedings  (I)  Dolgeshein  took  ,/|  = /2  and  ascribed  this  form  of  the  radi- 
ated pressure  to  Askarian.  It  is  called  the  Dolgeshein-Askarian  formula.  From  its  derivation  it 
is  seen  to  be  valid  only  when  kr^  > > 1 or  > > L/327t^.  Although  sin  H'/  M'is  a Faunhofer 
pattern,  the  factor  is  added  as  a Fesnel  (or  near-field)  correction. 

Heat  Model  IV  (Heated-Spot,  or  "Spike"  Model) 

Let  a particular  volume  f,)  in  infinite  space  suddenly  receive  a quantity  of  heat  which 
then  flows  outward  at  the  rate  //  = I,,  where  the  units  of  //are  calories  per  unit  volume 

per  unit  time  (N  • m/m-'  ■ s),  p is  the  density  of  the  medium  (N  • s/m'^),  is  the  specific 
heat  (m^  ■ s/K  ),  and  t/',  is  the  flux  of  temperature  (K/s).  The  equation  of  heat  conduction  to 
be  solved  is  then 

= D v^r  -f  V, 

rfl 


- — 


in  which  D is  the  diffusivity  constant  (units:  m ^/s ) and  T is  the  temperature  (K).  To  solve 
this  equation,  one  uses  an  appropriate  Green's  function  which  is  the  solution  of  the  auxili- 
ary equation 

1 


--b 


= - 4776  (r  -r„  ) ft  (/  ). 


" D di  ■’ 

Thus  the  solution  of  the  inhomogeneous  temperature  equation  is 


Now  (9) 


T =- f ill,,  f G/)(r,  r„|7.  7„)  H (t„.  i„)tn'„(x„)  . 

4 77  p C T) 

f//,  (/?,  r ) = I ^ ^ ^ • 
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Thus 


T = 


8ir 


VT 


l>Cp 


f'  ..f 


(Dt) 


3/2 


In  conformity  with  the  model  we  choose  Hlo  have  the  special  form  of  a point  source  in  shape 
and  time, 

H(ro.t„)  = £o8(r -r„)S(f -r„), 

in  which  the  units  of  Eg  are  calories,  i.e.  (N  m),  and  we  take  r^,  i„  to  be  at  the  origin.  Then 
the  history  of  the  temperature  distribution  at  space  point  r(  =|r|  ) and  time  t is 


Tir.  I)  = 


-rh%Di 


pC  (4nDr)^^^' 


We  now  turn  to  the  propagation  of  a thermoelastic  wave.  In  the  immediate  vicinity  of 
the  origin  the  temperature  rise  due  to  H causes  an  expansion  of  the  medium  in  the  amount  fiT 
(units  of  /i:m/m  K).  This  expansion  propagates  away  as  the  thermoelastic  potential  <1>  (units: 
m ^),  governed  by  the  equation 


V^(p  - 


1 

= -AttPT*.  T*  = T/4-n. 

/>2  ht2 


This  is  a small  amplitude  linear  equation  of  motion  in  which  c is  the  speed  of  propagation  of  the 
disturbance.  Since  the  Green’s  function  of  the  scalar  wave  equation  is  well  known,  the  solu- 
tion is  easily  derived  to  be 


1 r r h\Rjc -{t -ig)\ 

* - ^ J J >0 


dV., 


T(t„  I -R/c) 
£ 


or 


d>  (r,  t)  = 


PEp 

4npc„ 


f 


dV. 


g - R'^IADii  -RIC) 
" R[47tD(i 


■ E = I r -r„ 


The  integral  is  taken  over  a volume  of  radius  R =ct , with  center  at  r. 


As  it  stands  the  form  is  diiTicult  to  evaluate.  Nowacki  [10]  has  circumvented  these 
difficulties  by  solving  the  problem  in  a different  way.  He  begins  with  the  thermoelastic  equa- 
tion of  motion, 

/j.  V -1-  (\ +/U  )V  (divu) —■ yV  7 = pit 

(units  of  each  term,  N/m^),  in  which  u is  the  displacement  (m),  p.  is  the  second  Lame  con- 
stant for  an  isotropic  elastic  continuum  (units:  N/m^),  K is  the  first  Lame’ constant,  y is  the 
thermoelastic  constant  (N/Km^),  and  p is  the  medium  density.  He  then  seeks  solutions  of 
two  scalar  potentials  <l>,  and  ih  such  that 

u = Vd>  + curli//. 

Substitution  leads  to  an  equation  for  <l>: 

(K  + 2p)\7  (V^<P)  - yVT  = p-^  Vcl> 

dt^ 
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or 


where 


□ ?<1>  = mT. 


m2_r72  2 A +2// 

□ J _ V — and  m 

l>  A +2fjL 


= r_ 


Since  the  heat-conduction  equation  for  7"  is 


a ?T  = 


\ 2 =v2 


D 


D dt 


it  is  seen  that 


□ □ f (|>  = - 


Formally  this  can  be  written  as 


P 2 - O I 4>  = - ^ 


mU/ 

~D~ 

□ i 


Uf  = m(T  - S), 


where 


Hence 


T ^'1  ^ C ^'1 

/ = z r and  5 = — 


D aj 


J__d L_ii 

D di  f2  qP 


<l>  = mir  - S). 


Assuming  the  thermoelastic  continuum  was  free  of  all  strains  and  stress  for  t ^ 0,  one  can 
take  the  Laplace  transform  (written  with  an  overhead  bar)  and  obtain 


<1>  = m(T  - S). 


Since  □ — (1/c^)  and  03  = — (!/£))  B/Bf,  the  Laplace  transform  of  S' 

can  be  obtained  from  the  Laplace  transform  of  /"by  replacing  s/Din  the  latter  with 


f = C -R^/s,  0 
4nDR 


so  that 


Thus 


= -Q- 


AttDR 


-Ri/c 


. = 


+ - " 


(units:m  ’ ■ K ) . 


* = Q _ <’  -c  ‘ 

A-nR  c2 

v -^n 
C‘ 
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Laplace  inversion  leads  lo 


<1,  = !"  Q t'  (R  i)  + 0,  I < R/c 

AttR  ' /2  (R.i).i>  RIc  ■ 


Here,  J^  (R.i)  describes  the  heat  diffusional  effect  in  producing  displacement  potential,  and 
,/2  (R.i)  describes  the  wave  propagation  effect  (time  scale  much  faster  than  ./j).  Using  ./2  onl>, 
one  can  directly  show  that  the  displacement  potential  of  elastic  waves  produced  by  heat  is 


_ m Q 


-ii  -K  i ici  n _ I — 


The  pressure  wave  (shock  wave)  due  to  this  potential  is 


/) 

P = R r-  = p e 

a, 2 AttR 


Since  this  is  a transient,  it  will  be  convenient  to  form  first  the  pressure  impulse  (units: 
N • s/m ')  and  then  divide  by  some  average  time  Ar  to  obtain  the  time-average  pressure,  /> 
Thus 


•n  4ttR  O-  i 


^ ,,mQ  + C'^ 

AnRD 

Since  m = A/  (y  ),  and  A has  the  units  N/m  K,  m has  the  units  of  K that  is,  m is 
identical  to  fi,  the  coefficient  of  thermal  expansion.  Similarly,  since  one  can 

write 


J"  /i(i)  (Ji 


4nC],DR 


Thus  the  time-average  pre.ssure  will  be 


i-',.  ‘ ' 


AnC],  DR  Ar 


The  crucial  parameter  in  this  model  is  Ar,  the  duration  of  the  transient.  A discussion  of  Ar  is 
reserved  for  the  section  on  numerical  calculations. 

The  preceding  model  was  originally  derived  by  Bowen  |l),  and  is  called  the  Bowen  "hot- 
spot" model 

In  a later  paper  (ref.ll)  Bowen  reexamined  his  basic  formulation  in  the  light  of  contra- 
dicting experiments  on  received  acoustic  waveforms  which  displayed  "bipolar  transients"  (name- 
ly plus  and  minus  spikes)  rather  than  a mono  polar  transient  (or  plus  spike  only).  He  conclud- 
ed that  the  particular  solution  /^,  ( = <l> ) already  given  by  him  above  must  be  extended  to  a 
more  general  solution  containing  amii>lcnn’nian  so\uUor\s  (called  ih  ')  i.e. 


= 'l>^  (r.i)  -I- 
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in  which  i/<  ^ is  an  outgoing  wave,  and  i/i  ~ is  an  ingoing  wave.  The  solution  is  subject  to  the 
boundary  condition  that  at  the  origin  (r  ”0)  the  displacement  potential  •>•„„„/  must  vanish  for 
all  time  / ^ r/c,  viz. 


+ i/<  ■*■(0  - ci) 


+ ./(  " (0  + cr ) 


= 0 


in  which  d>^,  = (Kir)  In  addition  the  solution  must  satisfy  a set  of  initial  conditions. 

dd> 


=0;  ^ (^0) 


I i^rli 


= 0. 


By  use  of  this  formulation  Bowen  made  several  improvements  to  his  earlier  theory:  (I)  the 
infinity  at  r = 0 in  the  particular  solution  <l>  = d’  was  removed  (2)  the  discontinuity  di  r = tv 

' I 

in  <l>^,  (or  in  its  derivatives)  was  removed  (3)  the  radiated  pressure  impulse  /(/■,/)  ( = f p(r.  I') 


di'  = f>  dd>/9/)  became  l(r,i)  = 


£.,  H 

4 TT  dR 


[F(ei  - r.O  )1 


exp  [c(6)  + (a  - r)  ] 


erf 


. (ci  — r)  + 2c(t)  + i) 


[4D  ib  + /)]  ' 2 

in  which  h is  the  "initial  time"  (not  zero)  related  to  the  initial  size  < ^ of  the  heated 

sphere  by  the  formula  = 6Pb.  When  c(b/D)''^  > > 1 it  can  be  shown  that  the 

time  variation  of  Hr.t)  is  Gaussian, 

Toe  exp  [ — (cf  — r)^l4Db]. 

As  a result,  the  acoustic  pressure  is  bipolar,  i.e. 

_ 6/  iT  _ _ 2c(ci  - r) 
di 


/?  = OC 

dt 


ADb 


F. 


This  result  agrees  with  experiment.  The  first  moment  of  the  radiated  pulse,  namely 

J=  fit'  pirjldt  = r (i')  dt' 

bt 


can  be  shown  to  reduce  to  the  value 


y=  - 


47rCpr 


if  (again)  C(blD)^'^  >>  1.  This  agrees  with  the  magnitude  of  the  generalized  Westervelt- 
Larson-Hanish  model  explained  below. 

These  improvements  by  Bowen  move  his  model  in  the  direction  of  the  generalized  model 
used  in  the  calculations  noted  above.  Bowen  notes  again  that  the  pressure  itself  because  of  its 
ultrashort  time  duration  is  essentially  undetectable  directly  if  a single  hot-spot  is  in  question. 
However  7 itself  (as  defined  above)  is  finite,  and  certainly  the  superposition  of  many  hotspots 
is  detectable  if  £„  is  large  enough. 

Far-Field  Patterns  of  Pressure  Radiation 


It  was  shown  in  the  section  on  heat  model  HI  that  the  far-field  angular  distribution  of 
pressure  />(«)  is  proportional  to  the  Fourier  transform  of  the  spatial  distribution  of  heat  deposi- 
tion (designated  here  h(z))\  that  is. 
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I ig.  7 - Dcl'initions  of  aperture  function 
and  far  field  pressure  pattern 


p(u)  ~ J"  h(z)  e d:,  k = 2tt/K,  and  u = sin  W. 

The  angle  H is  defined  in  Fig.  7. 

Several  distributions  hi:)  may  be  pertinent  to  this  survey.  We  list  them  here,  together 
with  facts  on  sidelobes  of  the  developed  patterns.  In  each  case  we  take  the  length  of  the  depo- 
sition track  to  be  L and  choose  the  center  of  the  pattern  to  be  at  r = 0.  Patterns  other  than 
the  following  can  similarly  be  constructed  as  needed  to  conform  to  experimental  result,  or  to 
additional  modeling. 

Exponential  Heat  Deposition 

Let  /i(c)  be  proportional  to  e and  for  ease  of  computation  shift  the  line  c = 0 to  c = 

U2\  that  IS,  take  h(z)  = e r ^ 0.  Then,  as  sketched  in  Fig.  8 the  normalized  radiation 
pattern  is 

p(u)  . 

1 + / — sin« 

a V 

This  pattern  has  no  zero  crossing. 


l ig.  8 - Radiation  pattern  caused  by 
exponential  heat  deposition 
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Uniform  Heat  Deposition 

Let  hi:)  = reel  L,  having  unit  amplitude  between  —U2  and  +U2  and  being  zero  else- 
where. Then,  as  sketched  in  Fig.  9,  the  normalzied  radiation  pattern  is 

pi,)  ^ jin  U yx) 

TT  Lu/K 

The  first  zero  crossing  is  at  m = \/L,  and  the  first  sidelobe  is  at  ii  = 3\/2L  The  3 dB 
beamwidth  is  0.88  K/L,  and  the  first  sidelobe  is  —13.4  dB  down.  This  is  the  case  of  heat 
model  lil. 


t ig.  9 - Radiation  pattern  caused  b> 
uniform  heal  deposition 


Cosine  Heat  Deposition 

Let  hiz)  = cos  rect  Z.,  for  ^ c < with  hi:)  being  zero  elsewhere.  Then,  as 
sketched  in  Fig.  10,  the  normalized  radiation  pattern  is 

ptui 

- . TT  sin  IT  , sin  k|  Z'  Ny 

-T  — +-T“|’  \ 

where  


IT  = — — M - and  F = ii  + — 

K 2L  K 2L 


I'ig.  to  - Radiation  pattern  eaused  by 
cosine  heat  deposition 


The  first  zero  crossing  is  at  m = and  the  first  sidelobe  is  at  u and  has  a magnitude  of 
-23.5  dB. 


Gaussian  Heat  Deposition 

Let  hi:)  = e assume  no  limits  in  ;.  Then,  as 

sketched  in  Fig.  II,  the  normalized  pattern  is 


piu)  =<. 


I'ig.  1 1 - Radiation  pattern 
caused  by  gaussian  heat  depo- 
sition 


in  which  uy  is  its  standard  deviation,  namely  <r|  = K/ttL.  There  are  no  zero  crossings. 
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Parabolic  Heat  Deposition 

Let  h(:)  = |l  - (2.-//.)^]  reel  ~ J ^ ^ f 

Then,  as  sketched  in  Fig.  12,  the  normalized  beam  pattern  is 


pill)  r 


sin>l» 


— cos  >1) 


7X 


. >b  = 


kLii 


l ig.  J2  - Radiation  pattern  caused  by 
parabolic  heal  deposition 


The  first  sidelobe  occurs  at  ii  = and  its  magnitude  is  —22  dB.  The  3-dB  beamwidih  is 
1.16  K/L. 


Impulse  Response 

As  was  noted  earlier,  the  transient  acoustic  pressure  given  by 


pir.i)  =J  pir.oj)  ti((u)  e 


2it 

The  steady-state  pressure  is  a directional  function  of  the  spatial  coordinates.  A discussion  of 
spatial  dependence  has  just  been  presented.  We  are  interested  now  in  the  temporal  aspect  and 
omit  details  of  spatial  distribution  by  assuming  the  radial  pressure  to  be  spherically  symmetri- 
cal: 


pir.o))  =4 


, ikr 


To  find  A , we  shall  specify  some  radial  surface  velocity  of  a fictitious  sphere  (radius  a).  Let 
this  be  and  let  the  medium  of  propagation  be  homogeneous  and  isotropic.  Then 

lA 


kper^ 


1 — /Arl  k = m/c. 


and  on  the  surface  of  the  sphere 


V,.  )kpca^ 

Aia.o>) 


•lk(J 


/ (1  — ika) 


Hence,  the  acoustic  pressure  anywhere  is 


, , kpea^  e 

p(r.oj)  7~T~ 

/ ( 1 — ika ) 


/kir  ) 


V,.  (a.oj ). 


If  we  choose  a surface  radial  velocity  v^ia.i)  then 

\^(a.w)  = J c""'  til. 

Thus  the  transient  acoustic  pressure  radiated  into  the  far  field  is 


(/o> 

2n 
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As  an  example,  let  the  impulse  of  heat  deposition  be  a rectangle; 


v,(r)  = 
= 0 


0 < r < 7,„ 
elsewhere. 


The  Fourier  spectrum  is  known  to  be 


V,.  (a.iu)  = 


/ F., 


1 - c 


Uu  T,, 


By  direct  integration  one  obtains  the  transient  radiated  pressure  to  be 

p{r,i)  =0,  ; < /•  - a/c, 


= [/  - (/  -,/)  ( I 


< / < 


r —a 


= pcV„^ 


- (i,  </ ) If  - (r  -a),  (1  _ 


C 

t/)  U -(/  ■ 


+ T„. 


■I/)  i +7,1 


I > 


+ T.. 


Clearly  the  nature  of  the  response  will  depend  on  the  ratio  TJia/c).  Three  conditions  of  this 
ratio,  and  the  consequent  pressure  plots  are  shown  in  Fig.  13.  In  Fig.  13a  the  signal  p travels  a 
distance  cT,,  much  greater  than  the  radius  of  the  sphere;  both  positive  and  negative  impulses 
are  widely  separated  in  lime  and  are  therefore  equal.  In  Fig,  13b  the  signal  p travels  a distance 
<T„  which  is  approximately  the  same  as  radius  a.  The  positive  impulse  overrides  the  negative 
impulse,  reducing  the  latter.  In  Fig.  13c  the  travel  distance  cT„  is  smaller  than  a.  This 
corresponds  to  a spherical  radiation  suddenly  increasing  in  size  from  small  to  large.  The  nega- 
tive pulse  is  due  to  the  fluid  inertia  (moving  radially  outward)  at  the  time  of  cessation  of  ex- 
pansion, It  is  much  larger  than  the  positive  pulse,  because  the  larger  sphere  is  a more  effective 
radiator. 
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Another  time  impulse  is  the  triangle 

T 

V (t)  =.  nil,  0 < r < — . 

2 

= -nn  + 2 ^ < , < T„. 

= 0 : elsewhere. 

The  radiated  transient  appears  as  sketched  in  Fig.  14.  (Reference  [12]  gives  additional  de- 
tails.) 


I'ig.  14  - Impulse  response  caused  by  triangle  heat  deposition. 


Estimation  of  Time  Durations 


In  all  heat  models  the  estimation  of  time  durations  of  impulsive  heating  and  acoustic  ra- 
diation is  most  important.  Loosely  two  types  of  time-dependent  processes  must  be  considered: 
the  heat  flux  itself  (calories  per  second)  entering  the  energized  volume  of  the  liquid  and  the 
rale  at  which  the  thermoelaslically  induced  shock  waves  cross  the  heated  region  and  are  radiat- 
ed outward  as  sound.  The  two  lime  durations  are  explained  by  a simple  mathematical  model. 
It  was  noted  earlier  (for  heal  model  Ib)  that  if  the  heat  source  is  a delta  function  in  space  and 
lime  i.e.  ft  (/)  ft  (r) ),  the  acoustic  pressure  developed  is  a shock  transient: 


pir.t) 


1 

47r 


The  rale  of  energy  deposition  is  £„  fi(f)  or  £„/£,.  The  acoustic  time  history  is  S/dt  [£„  fi(r)] 
or  = £„/£]  7'2  estimates  of  Tj,  and  T^. 
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It  is  simplest  to  estimate  the  "acoustic  time"  characteristic  of  the  radiation  process.  Let  </ 
be  a characteristic  dimension  of  the  heated  volume  (e.g.  diameter  of  a cylinder  or  sphere).  By 
analogy  with  traveling  waves  in  a parallel-plate  duct  we  seek  a cutoff  frequency  of  such  waves, 
that  is,  the  frequency  associated  with  a mode  which  just  fails  to  propagate  down  the  duct.  If  / 
is  the  width  of  the  duct,  the  first  transverse  mode  (defined  as  reinforcement  of  incident  and 
rellected  waves  one  wavelength  K apart)  travels  at  an  angle  H such  that  21  cos  n = K (or  such 
that  A ^ 21).  Each  A is  associated  with  a particular  n.  The  largest  wavelength  for  the  first 
transverse  mode  is  A =2/at  W =0°,  which  means  the  wavefront  is  parallel  to  the  duct  wall  and 
forms  a standing  wave,  bouncing  back  and  forth  between  walls  at  one  wavelength  reinforce- 
ment and  not  propagating  down  the  duct.  The  cutoff  frequency  is  then  In  a 

similar  way  acoustic  reinforcement  in  a heated  volume  at  one  wavelength  separation  occurs 
when  A = 2d,  and  the  cutoff  frequency  is  = c/2(l.  The  choice  A = 2(1  is  equivalent  to  the 
wavenumber-diameter  product  kd  = rr  ( or  A =2vlK).  In  acoustic  theory  of  radiating  surfaces 
this  product  is  favorable  to  radiation  of  wavelength  A (favorable  in  that  the  normalized  radia- 
tion resistance  is  near  unity  and  the  reactance  is  low).  Alternatively  expressed,  it  is  favorable 
to  radiation  of  wave  numbers  near  ttIJ. 

The  cutoff  frequency  may  be  taken  as  the  bandwidth  A/ of  the  radiated  acoustic  power 
spectrum.  The  time  duration  of  the  associated  transient  is  At  = (If)  Thus,  roughly,  T2 
= ./„r*  • The  estimation  of  is  uncertain,  the  varying  values  put  forth  being  the  principal 
difference  between  models.  All  models  in  which  heating  of  the  fluid  is  the  source  of  sound  can 
be  reduced  to  the  schematic  formula 


= 


H E,, 


477  C.  7 


J_  _}_ 

t:  t: 


in  which  the  subscript  1 refers  to  heat  deposition  time  and  the  subscript  2 refers  to  acoustic 


transient  time.  We  assume  T2  =./2~  — .tm'  ■ Thus  T2  is  fixed  by  the  speed  of  sound  in  the 

liquid  and  the  characteristic  dimension  (=  r/)  of  the  model.  The  value  to  be  assigned  to  7j 
depends  on  the  model.  Four  models  and  their  assigned  values  are  as  follows: 


Heat  model  la  (infinite  with  space,  steady  state),  with 


\p\  = 


H En 


4rr  C„ 


277 


./2 


SO  that 


T,  = 


277 


where  is  unknown; 


Heat  model  lb  (infinite  space,  transient),  with 


IH  = 


4^  C',  I d 


il 


II  T2  = 77, 


so  that 


T,  — — where  T,  is  unknown 
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• Heat  model  III  (Dolgeshein-Askarian,  or  heated  rod,  model),  with 


IpI  = 


n Eo 


so  that 


4rr  kxl 


Tx  = 


2rr ./,  ./j 


^JTTTl 

2tt 


. /l  =/2. 


Heat  model  IV  (Bowen,  or  heated-spot,  model),  with 


= 


H Eo 


Air  C,  R 


DT, 


. 7;  = A T, 


so  that 


t; 


Thus  in  both  the  Dolgeshein-Askarian  model  and  the  Bowen  model  the  value  of  T,'  can  be  in- 
terpreted as  being  assigned.  For  the  other  models  we  can  estimate  in  the  following  (some- 
what arbitrary)  way:  We  assume  that  the  cosmic  particle  on  entering  the  liquid  is  decelerated  at 
a constant  rale  from  a maximum  velocity  v to  zero,  making  the  average  velocity  v/2.  The  time 
required  to  complete  the  transit  over  a length  / of  the  shower  is  therefore 


h 


2/ 

V 


For  a relativistic  particle  of  energy  E and  rest  mass  m^c}^,  Ci^  being  the  speed  of  light,  the 
velocity  is 


A typical  application  is  for  muons,  for  which  m„  = 1 x 10*  eV.  Since  cosmic  particles  of 
the  heavy  nuclei  type  have  energies  that  exceed  cosmic-ray  protons  (£  = 2.5  x lo’  eV),  the 
speed  of  muons  is  approximately  the  speed  of  light.  Thus 
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In  each  model  in  which  T\  is  not  assigned  we  choose  some  characteristic  length  / to  estimate 
T]  by  this  formula. 


The  arbitrariness  of  this  estimate  is  apparent.  From  tentative  experiments  [1]  and  other 
models  [7]  it  seems  more  plausible  to  lake  to  be  about  l//s,  which  is  some  2 or  3 orders  of 
magnitude  larger  than  the  values  predicted  above.  This  estimate  is  sharply  different  from  what 
IS  predicted  from  bubble  models,  from  which,  as  discussed  latter  in  this  report,  T,  appears  to 
be  of  the  order  10  ~".s. 


Numerical  Calculations 


The  models  listed  above  will  now  be  numerically  evaluated.  In  each  case  we  will  present 
•he  absolute  magnitude  of  the  pressure  field  in  the  direction  in  which  it  is  a maximum,  at  a 
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nominal  distance  ol  I ni  from  Ihe  origin  (or  radiation  center),  and  per  unit  energy  input  (in 
electron-volts).  When,  as  discussed  in  the  preceding  subsection,  the  deposition  time  is  not  as- 
signed, we  will  assume  the  interaction  shower  to  be  a cylinder  a cm  in  radius  and  /.  m long, 
the  values  ol  a and  /.  being  assumed  in  particular  cases. 

Model  la 

In  model  la 


l/'l  = 


^7T  C. 


VT” 


VI  -f  cos  ^ 0 

We  take  /„.S’  = t and  = 2tt  where  /i  = , and  we  calculate  the  maximum 

at  broadside,  which  is  at  0 =90°,  To  calculate  T|,  we  assume  the  shower  length  L to  be  10  m 
(ll.  Thus 


21. 


2 X 10  m 


T, 

‘ / 0.3  X 10*^  m/s 


= 67  IIS. 


To  estimate  Zj,  we  will  take  /i  — i,2(/,  where  </  is  the  diameter  of  the  cylinder.  The  diameter 
depends  on  the  energy  and  is  esiimateil  to  be  between  1 and  12  cm.  Let  us  assume  (/  = 3 cm, 
so  that 

1.5  X lO^'cm/s 


/i  = 


3 cm 


= 25  A//.- 


and 


Thus  at  1 m 

Izi  =__/^ 

£■ 


7,  =/2~'  MOms. 


47r  C;,  T,  Tj  {r| 

1.4  X 10  K X 1.6  X 10  y/cl 

471-  X 4.18  X 10-*  ( m^/s^  K ) x 40  x 10  s x 1 m x 67  x 10  '''  s 
= 1 .6  X 10  '^N/m^-  eV  = 1 .6  x 10  dynes/cm  ^ • eV'  . 

Hence  for  Model  la 


~ Pa/eV  rms  at  I m from  the  axis  of  a rod  10  m long  and  3 cm  in  di- 

ameter, for  steady-state  conditions  at  25  kHz,  with  7,  =67  ns. 


The  total  power  radiated  in  the  case  kl  > > I is 


(/„  .V)‘ 


acouslic 


16  II  (■; 


Setting  7„.V  = a ' = 7.,  we  estimate 


(/<  tmstn 
)2 


= 2:r  X 2.5  x IQ-*  x (|  4 c ;0  4)7  ^ ^ .,5 

7.  !6p(;  10  X 16  X 10^  x (4.18  x K)’)^ 
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Hence  lor  model  la  the  loial  radiated  power  is 


*^atousiic  ~ '-'0  ^ ffis  from  a rod  10  meters  long  at  25  kHz. 


Model  Ib  (One-Shot  Siren) 


For  model  Ib  the  time  history  of  the  transient  radiation  is 


p(r,i)  = 


li  11  o (/„  S) 


4n-  C + A'^cos^«)|rl 


(>  cos  11 


— A'  cos  H sin  11 


,-iiL 


n 


(I  - R/c)  - T/2 


. K = 11/f. 


At  broadside  « =90°;  Hence  the  maximum  amplitude  is 

11 


IpI  = 


477  c„  H T, 


As  before  we  take  11  Fj  = tt,  so  that 


IH  = 


4 C,  IH 


_1_  J_ 

Fi  T, 


This  is  77  times  larger  than  model  la;  that  is,  for  model  Ib 


/>n,ax  = 5.03  X 10  Pa/eV  at  1 m from  the  axis  of  a rod  10  m long  and  3 cm  in  diameter, 
for  transient  conditions  at  center  frequency  of  25  kHz,  with  T\  = 67  ns. 


Model  HI 

In  the  Dolgeshein-Askarian  model 

, , h t-o  1 1 sin  W" 

P (r  7))  = =r — 

^ 2 C , VZT;  if 

On  the  axis  at  broadside  sin  WIW  = \.  Taking  ,/]  =/2,and  F,  = F2,  one  arrives  at 

l/^l 

2 C yj  Ll\ 

In  the  near  field,  at  i\  = 1 m, 

IH  1.4  X 10  (K)  X 1.6  X 10  -I'^J/eV 

F„  2 X 4.18  X 10^  (mV.s^  K)  (10  x D'^m  (4  x lO'^)!^!- 
Hence  for  model  II! 


F2F2 


1/>|/F„  = 5.2  X 10  Pa/eV  rms  at  1 m on  the  axis  of  a rod  10  m long  for  steady- 
state  conditions  at  25  kHz. 


■j 


» 

( 

[ 


36 


/ 


NRL  REPORT  8150 


Model  IV 


In  the  Bowen  model  Ihe  average  pressure  radiated  is 


47rC;,/)/?AT 


Thus 


> a\ 

~e7 


T,  = 


D 


The  effective  deposit  time  is  given  by 

^ D _ 1.43  X 10  m/s‘ 


and  7'2  = At. 


1 


(1.5  X 10^)2  m/s^ 


= 6.35  X 10 


-'‘’s. 


and  the  interval  between  compression  and  rarefaction  is 

h 


= 4 X 10  “5  s. 


J 


Hence 


1.4  X 10  K X 1.6  X 10  J/eV 


47r  X 4.18  X 10^  (m^/s^  K)  x 6.35  x 10  s x (4  x 10  "5)  s x 1 „i 


= 1.67  X 10 


-to 


N/m^  ■ eV 


Thus  for  the  Bowen  model 


pJE„  =1.67  X 10 
corresponding  to  25  kHz. 


-10 


Pa/eV  at  1 m from  the  origin  of  a spike,  averaged  over  AO/xs. 


Ambient  Noise  and  Molecular  Agitation  Noise 


Below  35  kHz  the  noise  in  the  ocean  at  the  lowest  measured  level  is  the  ambient  noise, 
in  distinction  to  molecular  agitation  noise.  The  lowest  limit  of  prevailing  (ambient)  noise  has 
been  compiled  for  all  the  oceans  by  Wenz  (ret.  13).  Table  1 is  a brief  listing  of  sound-pressure 
spectrum  level  in  a 1 Hz  band  in  various  units  taken  from  this  reference.  Above  35  kHz  the 
noise  L,  (units:  dB)  due  to  molecular  agitation  is  predominant.  It  can  be  calculated  from  an 
empirical  formula  at  frequency  ./ (units:  Hz), 


75  + 20  log|o../,  dB  re  1 ^Pa 


(14).  At  25  kHz 


L,  = -75  -I-  20  log|„  25,000  = 13  dB  re  1 ^Pa.  = -87dB  re  //bar  . 


This  is  5 dB  less  than  ambient  (Table  1).  In  making  our  calculations  we  will  use  the  (rough) 
noise  figure  of  7.9  x 10  Pa  in  a 1 Hz  band. 


To  process  the  acoustic  transients  in  their  entirety,  one  must  arrange  to  have  sufficient 
bandwidth  in  the  receiver  circuit.  We  noted  earlier  that  the  smallest  time  intervals  to  be 
recorded  are  ol  the  order  of  40  //s.  Roughly  we  take  the  required  bandwidth  to  be  10“*  Hz 


Hi 


/ 1 
i 

A 
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Table  I - Ambient-Noise  Sound  Pressure 


Spectrum  Level  in  a I -Hz  Band 


f 

(Hz) 

tx  Pa 

dU  re 

20 /a  Pa 

dB  re 

1 n bar 

1 

28,300 

63 

-1 1 

10 

1,120 

35 

-40 

100 

100.2 

14 

-60 

1,000 

22.4 

1 

-73 

10,000 

10 

- 6 

-80 

25,000 

7.y 

- 8 

-82 

30,000 

6.3 

-10 

-84 

Hence  the  total  noise  spectrum  level  against  which  the  signal  spectrum  level  to  be  detected  is 
(again  roughly)  7.9  x 10  x VlC*  = 7,9  x 10  Pa  at  25  kHz  (say  approximately  10 
Pa).  Here  the  spectrum  level  is  a.isumed  flat  over  the  band.  However,  if  we  do  not  desire  to 
know  the  entire  signature  of  the  transient  radiation,  the  receiver  bandwidth  can  be  much 
smaller.  The  theory  of  signal  detection  in  noise  allows  the  choice  of  receiver  bandwidth  to  op- 
timize the  probability  of  detection  against  the  probability  of  false  alarms  relative  to  a threshold. 
Thus  the  bandwidth  is  fixed  only  to  the  extent  that  probabilities  and  thresholds  are  selected 
subject  to  the  goals  of  the  signal  processing  itself.  We  continue  the  analysis  on  the  arbitrary- 
basis  of  a 1-Hz  band  and  a flat  \0‘’-Hz  band. 

Limits  of  Detection  Range 


Local  noise  (ambient  or  molecular)  limits  the  detection  of  a distant  signal  to  certain 
ranges.  We  choose  a criterion  that  a signal  is  detectable  if  its  absolute  level  is  equal  to  or  is 
above  the  sound  pressure  spectrum  level  of  local  noise  over  the  bandwidth  of  the  receiver.  In 
practice  the  threshold  is  set  higher  than  the  local  noise  level,  but  here  we  are  seeking  orders  of 
magnitude.  Applying  our  criterion,  we  see  that  the  limit  range  of  detection  (/?/,„,„)  in  units  of 
meters  per  eV  may  be  obtained  from  the  previous  numerical  calculations  by  use  of  the  formula 


R 


n 

liniif 


X 


fifr'  ^ = 


\pl ^o\  K =1 


in  which  //depends  on  the  model  and  SIT' is  the  bandwidth  of  the  receiver  in  hertz  (relative  to 
power).  We  choose  the  local  noise  to  be  in  a 1-Hz  band  at  25  kHz:  7.9  x 10  Pa/ (Hz) ' 
Applying  the  formula  to  selected  models  we  have  the  following.  In  these  calculations  we  as- 
sume attentuation  is  due  only  to  spherical  spreading  and  neglect  absorption. 
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Model  la 


In  model  la,  n = I and  \ plE„\ ^ =1.6  x 10  Pa/eV.  Thus 

x 2 1-6  X 10  Pa/eV  x ] n 

7.9  X 10  Pa/Hz'  ^ 
-II  p-  'Hz)'  ^ 

“ r 


' lima 


= 2.03  X 10 


from  which  we  obtain 


R 


liniif 


= 2.03  X 10 


I ni 
eV 


for  flH'  = 1 Hz 


and 


= 2.03  X 10  for  BW  = lO^*  Hz. 

e V 


Model  fb 


In  this  model  lb  ii  = 1 and  the  maximum  amplitude  of  transient  pressure  measured  at  25 
kHz  is  5.03  x 10  **  Pa/eV.  The  range  limit  is  therefore 


I 2 _ 5.03  X 10  -'<>  Pa/eV  x i m 


= 6.37  X 10 


7.9  X 10  Pa/Hz  ' ^ 
, 1 2 


-I I mHz 


eV 


from  which  wc  obtain 

=6.37  X 10  m/eV  for  SB  = 1 Hz 


and 


'limit 


= 6.37  X 10  m/eV  for  BW  = lO^*  Hz, 


Model  III 

In  model  111, /;  =yand|/7/fJ  =5.2  x 10  Pa/eV  x m'2  The  range  limit  is 
J 2 „ Dij,i  2 5.2  X 10  “’’Pa.'eV  x m'  ^ 

\limii  X 7 — 

7.9  X 10  Pa/Hz  ' ^ 


= 6.58  X 10  -'“m'  2 . Hz'  2/eV. 


trom  which,  lor  points  of  observation  in  the  near  field, 

'I'/.i,  = 6.58  X 10  -'•»  for  BW  = 1 Hz 
eV 


and 


= 6.58  X 10  for  SB'  = lO'*  Hz. 

eV 


In  the  far  (icld  replaces  ^ . 
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Model  IV 

In  model  IV,  /;  = I,  and  = 1.67  x 10  Pa/eV  x 1 ni.  The  range  limit  is 

Ri  X SIV''7  = 1-67  X 10  dynes/cni^/eV  x I m 


7.9  X 10  ■<’  Pa/Hz'  ^ 
= 2.11  X 10  m • Hz'  ^/eV. 

From  '*(hich 

^hmn  — 2.1 1 X 10  m/eV  for  BW  = 1 11/ 


and 


R,„„„  =2.11x10  ’ m/eV  for  BW  = 10  sup  4 Hz. 


MICROBLBBLE  MODELS 

The  mechanism  of  the  generation  of  sound  by  high-energy  particles  in  water  described  in 
the  preceding  section  is  that  of  a thermal  shock  of  an  elastic  medium.  We  consider  next  a 
different  mechanism,  in  which  sound  is  produced  by  formation  and  collapse  of  microbubbles. 
The  theory  of  this  process  originated  with  Selle  llSl,  who  hypothesized  that  cavitation  nuclei 
in  liquids  may  be  formed  in  part  by  irradiation  by  cosmic  particles.  In  several  subsequent  arti- 
cles Sette  and  Wanderlingh  [16,17]  elaborated  on  this  hypothesis  both  with  theory  and  experi- 
ment and  applied  it  to  explain  the  phenomena  of  bubble  chambers.  Although  they  did  not 
discuss  the  connection  between  microbubble  formation  and  collapse  to  acoustic  noise  genera- 
tion, it  is  a direct  matter  to  relate  these  two.  We  will  make  this  connection  later,  but  first  we 
will  review  the  theory  of  bubble  nucleation  to  lay  the  basis  for  numerical  calculation. 

Bubble  Nucleation 

Historically  the  first  significant  paper  on  bubble  nucleation  by  high-energy  particles  in 
liquids  was  that  of  Seilz  [181  in  connection  with  the  theory  of  bubble  chambers.  He  hy- 
pothesized that  most  of  the  bubbles  are  nucleated  by  energetic  free  electrons  resulting  from 
collisions  of  the  particles  with  the  molecules  of  the  liquid  medium  in  the  bubble  chamber. 
These  electrons  are  decelerated  with  extreme  rapidity,  producing  localized  hot  regions  or  "ther- 
mal spikes"  which  explode  into  bubbles  of  larger  than  critical  size  (that  is,  the  minimum-size 
bubble  is  not  in  equilibrium  with  its  surroundings,  which  permits  growth  by  vaporization  from 
Its  walls  into  in  cavity  at  the  temperature  of  the  cavity).  The  time  of  bubble  creation  by  explo- 
sion he  calculated  to  be  of  the  order  of  O.l  to  O.Ol  ns.  The  bubble,  after  being  created,  is 
thought  to  expand  by  withdrawing  heat  from  the  immediate  surrounding  Iluid.  The  maximum 
bubble  size  attained  is  limited  by  the  initial  energy  of  explosion,  by  heat  diffusion  away  from 
the  bubble,  and  by  depletion  of  local  permanent  gas  (air)  in  equilibrium  with  the  surrounding 
medium. 

The  energy  for  growth  o(  bubbles  in  the  the  bubble  chamber  is  wholly  supplied  by  the 
medium.  In  contrast  to  the  theory  of  bubble  chambers  the  production  of  microbubbles  in  or- 
dinary liquids  by  high-energy  particles  requires  that  both  creation  and  growth  be  energized  by 
the  kinetic  energy  of  the  particle,  no  energy  is  supplied  by  the  medium.  In  bubble  chambers 
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the  temperature  of  the  medium  is  locally  near  the  boiling  point.  When  a high-energy  panicle 
ionizes  the  medium,  the  particle  supplies  energy  to  create  an  embryo  of  critical  size.  Growth 
then  proceeds  by  evaporation  of  the  medium  at  the  bubble  wall,  the  heat  being  supplied  by  the 
medium  to  the  embryo  in  an  isothermal  process.  In  contrast,  when  a high-energy  particle 
enters  an  ordinary  liquid  whose  temperature  is  well  below  the  boiling  point,  there  exist  no  em- 
bryos with  critical  radius.  First  a limited  region  of  medium  must  be  heated  above  the  boiling 
point,  and  in  it  an  embryo  must  be  formed  which  has  the  critical  size  as  determined  by  the  ini- 
tial temperature  of  explosion.  Then,  as  the  embryo  expands,  it  draws  heat  from  its  heated  sur- 
roundings, cooling  down  a local  heat  pocket  ("spike")  in  the  medium.  As  the  temperature  of 
the  heated  pocket  falls,  the  theoretical  critical  size  of  a microbubble  increases  (Fig.  15)  at  the 
same  time  as  the  radius  of  the  actual  expanding  bubble  (of  initial  energy  Ej  increases,  the 
latter  being  larger,  but  the  difference  between  them  eventually  growing  smaller  because  of  ex- 
haustion of  the  original  heat  of  deposition  by  heat  transfer  and  gaseous  diffusion  At  some 
lower  temperature  the  bubble  radius  becomes  critical,  then  growth  ceases  (point  A).  If  £„  is 
large  enough,  the  critical  radius  is  not  reached  at  any  lower  temperature:  the  bubble  keeps  on 
expanding  (curve  £„)  until  stopped  by  exhaustion  of  the  initial  energy. 


t ig.  15  - Theoretical  critical  radius  vs  tcmperatuie  of  medium 
(I  ig.  7 of  Ref.  16) 


A more  intimate  description  of  bubble  nucleation  has  been  developed  by  Sette  and  Wan- 
derlingh  [11).  From  various  experiments  and  thermodynamics  they  have  constructed  the  fol- 
lowing theory  of  the  energy  balance  in  bubble  formation.  The  generation  of  a bubble  of  ra- 
dius R requires  an  energy  E(R)  which  is  the  sum  of  two  terms:  the  energy  required  to  form 
an  interface  of  given  surface  tension  fixed  by  the  temperlure  of  the  bubble  and  the  energ\  re- 
quired to  evaporate  water  from  the  interface  into  the  bubble.  Thus 

E(R)  = AttR^.t*  + iti  A/y. 

in  which  ir  * is  the  surface  energy  density  m is  the  mass  of  liquid  (water)  evaporated  in 

UK 

grams  and  \H  is  the  heat  of  vaporization  in  J/g.  A calculation  based  on  this  equation  shows 
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lhai  hundreds  of  millions  of  electron-volts  (1  eV  = 1.6x  10  are  needed  to  produce 

bubbles  of  raduis  10  cm  and  larger  and  R is  the  radius  in  m.  Experiments  show  that  only 
10  MeV  (approximately)  are  needed,  thus  contradicting  the  calculation.  The  discrepacy  is  ex- 
plained by  noting  that  the  role  of  dissolved  gases  has  been  neglected.  Actually  the  heat  A//  is 
the  sum  of  the  heat  associated  with  the  transition  of  dissolved  gases  from  the  liquid  into  the 
bubble  cavity  and  the  heat  of  evaporation  of  water  from  the  bubble  wall.  The  first  heat  is 
negative  (energy  is  released  when  gaseous  air  comes  out  of  solution),  and  the  second  heat  is 
positive.  fiiR)  is  always  positive  and  goes  to  zero  at  the  critical  point  of  the  liquid.  Below  450 
K the  gas  diflusion  predominates  and  AW  is  negative  [Fig.  5 of  Ref.  16].  The  importance  of 
this  case  requires  additional  comment.  If  one  rewrites  the  energy  balance  in  the  form 


£'(«) 


-I-  A//  (RP„  + 2<t)/3RT 


in  which  the  braces  include  the  total  surface  energy,  AW  is  the  sum  of  gas  and  vapor  heats, 
and  iR/’,,  + 2<r)l3RT  is  the  mass  in  grams,  assuming  a perfect  gas  (R  being  the  gas  con- 
stant), and  if  this  equation  is  plotted  as  the  expression  in  braces  vs  absolute  temperature,  it  will 

be  seen  [Fig.  6 of  Ref.  16]  that  the  unit  surface  energy  (the  entity  in  braces)  required  to  form  a 
bubble  of  radius  R = lOO/u//;  is  negative  below  410  K and  positive  above  it.  This  means  that 
below  410  K the  heat  generated  by  the  freeing  of  permanent  gases  from  solution  into  the  cavi- 
ty is  more  than  adequate  to  evaporate  the  medium  at  the  wall  and  sustain  surface  tension  of  a 
lOO/j/n  bubble  and  above  410  K it  is  inadequate  and  energy  must  be  supplied  from  the  cosmic 
particle.  Thus  at  600  K,  wnich  is  the  threshold  for  the  sure  formation  of  a cavitation  nucleus, 
an  external  source  of  energy  will  be  required  to  form  a 100/a/?;  bubble.  In  Ref.  this  energy  is 
estimated  to  be  4 to  5 MeV  for  a 70/a/;;  bubble,  both  as  to  its  initiation  and  its  "indefinite" 
growth.  Energies  less  than  this  create  bubbles  of  lesser  size  whose  growth  ceases  when  a criti- 
cal radius  is  reached.  This  phenomenon  of  bubble  growth  and  cessation  is  vital  enough  to  re- 
quire additional  comment.  It  can  be  pictured  as  follows  (Fig.  16).  A marble  is  pitched  up  the 

slope  with  an  initial  energy  E, , < 4 MeV.  This  brings  the  marble  to  radius  /?|,  where  motion 

ceases  (the  bubble  stops  growing).  If  the  energy  is  E2,  also  less  than  4 MeV,  the  marble  stops 
at  R2-  However,  if  the  energy  is  4 to  5 MeV  the  marble  rolls  over  the  threshold  and  continues 
on  "indefinitely"  (the  bubble  expands  indefinitely  by  absorption  of  heat,  if  available  from  its 
surroundings). 


I he  phenomenon  of  zero  or  negative  energy  required  to  form  bubbles  at  particular  tem- 
peratures may  lead  to  the  conclusion  that  there  is  no  limit  to  the  number  of  bubbles  that  can 
be  formed.  However  each  bubble  formation  is  based  on  the  availability  of  an  embryo  or  nu- 
cleus. The  probability  of  finding  such  a nucleus  depends  strictly  on  the  initial  temperature  of 
the  heated  region.  Above  600  K this  probability  is  near  unity  (assuming  atmospheric  pres- 
sure), and  below  600  K the  probability  becomes  negligible  [17],  Hence  there  must  be  enough 
energy  to  bring  a sufficient  volume  of  water  to  600  K.  Making  the  assumption  that  the  water 
so  heated  is  saturated  leads  to  an  initial  water  pressure  (in  a local  area)  of  12.2  MPa  (1786.6 
psi)  abs  (from  steam  tables).  After  expansion  to  ambient  pressure  (or  below)  the  bubble 
reaches  a final  radius  /?niax  value  is  determined  by  the  initial  energy  of  the  interaction 

products. 

The  expansion  to  lower  pressure  occurs  so  rapidly  it  can  be  thought  of  as  an  explosion. 
It  is  pictured  in  this  way:  Water  at  room  temperature  is  heated  at  nearly  constant  volume  to 
600  K (620°F).  as  if  the  water  were  in  a tiny  steel  tank.  The  momentary  existence  of  this  tank 
IS  attributed  to  the  high  rate  of  deposition  of  the  energy  from  the  cosmic  particle  (or  its  in- 
teraction products).  When  sufficient  heat  is  added  to  bring  this  liny  volume  of  water  to  satura- 
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lion,  the  saturaiion  pressure  in  the  fictitious  tank  is  12.2  MPa,  as  noted  above.  At  the  instant 
of  completion  of  deposition  the  lank  explodes  into  a bubble  which  expands  until  the  energy  of 
deposition  originally  supplied  is  exhausted.  The  expansion  is  nearly  isothermal  as  the  bubble 
draws  energy  from  the  heated  tank.  Toward  the  end  of  the  cycle,  when  no  more  heat  is  avail- 
able, the  bubble  can  continue  to  expand  briefly  in  an  adiabatic  process  down  to  ambient  pres- 
sure or  below,  the  driving  force  being  the  residual  inertia  of  the  water  immediately  surround- 
ing the  bubble. 


Theory  of  Forster  and  Zuber 

The  preceding  description  of  bubble  formation  and  growth  has  yielded  two  important 
parameters;  the  driving  initial  temperature  and  the  final  (or  maximum)  radius.  The  driving 
temperature  (here  taken  to  be  600  K)  initiates  the  explosion,  and  the  initial  input  energy 
determines  the  final  radius  of  the  bubble.  The  lime  history  R{t)  of  the  radial  expansion  is 
however  not  simple.  Numerous  authors  have  investigated  this  knotty  problem.  We  will 
proceed  here  to  sketch  an  approximate  theory  based  on  the  work  of  Forster  and  Zuher  |19l. 
which  has  been  selected  because  of  the  clarity  with  which  it  brings  the  intractable  parts  of  the 
theory  into  focus. 
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The  equation  of  radial  motion  of  a bubble  wall  begins  with  an  extended  version  of  the 
early  formulation  of  Rayleigh.  To  set  this  up  we  note  that  the  rate  of  change  of  volume  veloci- 
ty of  a bubble  of  radius  R is 


~-^nR^  =4ir  2RR^+R^r]. 

ill  ill  r t 


di\dt  2 \ ' ' 

If  one  assumes  the  liquid  is  incompressible  and  then  applies  Bernoulli’s  law  of  the  flow  of  an 
incompressible  fluid,  one  arrives  at  an  equation  of  motion 

R'^R  + 2RR^  R*R^  P,n,  - Pcxi 


in  which  r is  any  radius  '^R  and  the  right-hand  side  is  the  difference  between  the  internal 
pressure  and  the  external  pressure  (divided  by  p^).  Setting  r = R and  writing  the  internal 
pressure  as  the  sum  of  the  vapor  pressure  p,  and  the  permanent  gas  pressure  and  the  exter- 
nal pressure  as  the  sum  of  the  surface-tension  pressure  p„  and  the  hydrostatic  pressure  p^, 
one  can  then  write 


RR  + ^ r'^ 


Pv  + /’v  - Poo  + • 


To  simplify  matters,  one  first  neglects  the  internal  gas  pressure  and  then  relates  the  quantity 
■^P  ” /’v  “ Poo  to  tfie  change  in  temperature  AT  by  means  of  the  Clausuis-Clapeyron  equation 

T(v^  - V,) 

in  which  and  denote  the  specific  volume  of  the  vapor  and  fluid  respectively,  L is  the  la- 
tent heat  of  vaporization,  and  AT  is  the  rise  in  temperature  above  T.  The  temperature 
difference  AT  is  determined  by  the  solution  of  a problem  involving  heat  conduction  across  a 
moving  surface  of  evaporation  (in  spherical  geometry).  This  solution  may  be  found  in  Ref. 
[15].  If  the  heat  input  to  the  fluid  is  ^calories  per  unit  volume  per  unit  time,  then  for  thermal 
dififusivity  D 

K -rl  ^ H 1 1 1 \ i 


AT(/)  = 


C/inD)' 


Jit)  -f  At  + 


Here  and  p,  are  the  densities  of  vapor  and  liquid  respectively,  D is  the  thermal  diflfusivity  of 
the  fluid,  Jir.r.i)  is  an  integral  functional  arising  from  the  solution  of  the  heat-conduction 
equation  noted  in  Ref.  (201,  At  is  the  superheat,  and  C,  is  the  specific  heat  of  the  liquid  cal/gm 
K.  Subsitituting  all  these  expressions  into  the  equation  of  motion,  one  obtains 


RR  + ^ R^ 


- li  + yJil)  = Q*i, 


where 


r T(V„  - V,  ) 


p/T(v^  -v,)(nD) 


fijC/Tiv^  V/) 

This  is  an  integral-differential  equation  in  Rii):  Jit)  is  a functional  of  R{r).  A more  reveal- 
ing equation  can  be  obtained  by  introducing  the  concept  of  critical  radius.  This  is  done  by  not- 
ing that  for  superheat  At  there  is  a corresponding  pressure  rise 
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^Ps  = 


L\t 


T(v^  - Vj) 


= lit'/- 


The  bubble  radius  corresponding  to  a surface  tension  ir  at  pressure  is  called  the  critical  ra- 
dius 

p — = 2(r  _ 

" Ap,  lip,  fi  ■ 

Dividing  the  equation  of  motion  by  more  useful  form, 


+ — 1^  - — + BJ  = qi. 


where 


A = 


2(t 


P/^cr, 


. B = 


R‘ 


.^.r  = 


R 


Ktn, 


R. 


crn  *'crii 

Because  of  the  functional  integral  Jit)  the  solution  of  this  equation  is  intractable  as  it  stands. 
Hence  a general  formula  r{i)  for  all  time  i is  not  available.  Forster  and  Zuber  demonstrate 
that  for  very  small  bubbles  one  can  neglect  the  water  inertia  during  growth.  They  then  reduce 
the  problem  to  the  solution  of  an  equation  of  the  Volterra  type, 

r - 1 , - , 1 f'  r'(t)r(t')  A 

+ </(')  = TTT  I 
C r •'  0 


(/  - /') 


dt'.  c*  = 


in  which  Jh)  has  been  reduced  by  various  arguments  to  the  integral  shown.  Even  this  form  is 
intractable.  However  an  upper  limit  for  r(t)  may  be  obtained  by  using  the  mean-value 
theorem  to  calculate  the  integral  (which  then  becomes  Using  this  approximation  and 

assuming  qU)  is  small  relative  to  the  other  terms,  one  finally  arrives  at  a solution  for  rii)  im- 
plicitly given  by 

. ..  J2 

r - 1 


r + In- 


- 1 


= C*f ' 


I > 


The  initial  condition  is  r = r,  for  / = r^/C"^.  To  start  the  numerical  calculation,  one  assumes 
r.  is  slightly  greater  than  l,say  r,  =1.01.  Since 


2,r  pfTiv^  - v,){nD)^%R^, 


C*  = 


P l^cril 


L^P,,R,r, 


_ ,,,\T{nD) 

^irii  K 

we  can  interpret  AT  to  be  the  initial  superheat.  Hence 


C*  = 


\tC,AttD) 

R^.^i,Lv, 


1.2, 


This  is  a useful  formula  for  estimating  the  numerical  value  of  the  right-hand  side  of  the  solu- 
tion. To  obtain  a feeling  for  the  magnitude  of  the  terms  and  parameters  involved  we  will  calcu- 
late two  cases  as  follows. 

Case  I.  The  liquid  (water)  is  superheated  5°C  above  boiling  at  atmospheric  pressure.  We 
desire  to  find  and  the  time  required  to  form  a bubble  of  radius  R = 20  pm. 
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Solution  to  Case  J:  The  thermodynamic  chart  for  water  shows  that  at  T = 373  + 5 = 
378°K,  - 27  ft^/lb,  v,  - 0.02  ft^/lb,  L » 974  Btu/lb,  C,  = 4.18  J/g/  K,  and  D = 1 43 

X 10  cm  Vs.  Then 

(974  X 2.32)  (J/g)  X 5 “A-  X 0.02  x 62.4  cm  Vg 


LStv. 


- V,) 


378  “A:  I (27  - 0.02)  x 62 .4 1cm  Vg 
= 2.3  X lO^cmVs^ 


and 


« =—  =2frv/  = 2 X 50  X 0.02  x 62.4  = 124.8cm  Vs 


Therefore 


R - SL 

^cnt 


124.8cm  Vs 


3/.  2 


^ 2.3  X lO^cmVs^ 


5.6  X 10  '*  cm  = 5.6  nm. 


From  this 


St  (ttD) 


(■  STCfrr  X 1.43  X 10  ~^)‘/Vm/s'^^  x 4.18  (J/g  K)  x 27  x 62.4cm  Ve 
5.6  X 10  -^cm  X 974  x 2.32  (J/g)  x 0.02  x 6.24  cm  Vg 
- 1494  s 


Since  r — = 20  fjim/5.6  iim  = 3.57,  we  find 


r + In-^^ \ C*t'^^ 


or 


3.57  + In  = 1494/ 

or 

/ = 37  /US. 

This  is  approximately  the  result  that  Forster  and  Zuber  show  in  their  Fig.  1,  where 
/?  - 2 X 10  ^ cm  corresponds  to  a formation  time  of  about  30/us.  Note  that  / depends  on  C * 
inversely;  that  is,  the  higher  the  superheat,  the  faster  the  time  of  formation  of  a bubble  of 
specified  3ize. 

Case  2.  We  now  calculate  the  case  of  a cosmic  particle  (or  interaction  product)  with 
enough  energy  to  form  a 600-K  spike.  We  desire  to  find  the  time  required  to  expand  an  em- 
bryo to  a radius  of  70  fim. 

Solution  to  Case  2.  The  superheat  is  227‘’C.  From  steam  tables,  using  vapor  at  15  psi  and 
620“F,  and  fiuid  at  15  psi,  212‘“F,  we  have  - 43.6  ft^/lb  - 2.72  x lO^cmVg,  v,  = 0.02 
ft  /lb  - 1.25  cmVg,  latent  heat  of  vaporization  is  1345  - 180  - 1165  Btu/lb  = 2.7  x 10 
erg/g.  The  specific  heat  at  constant  pressure  is  4.18  x 10 ’ erg/g  K,  and  the  diffusivity  is  I 43 
X 10  "'J  cm  Vs.  Then 
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fi  = 


L\rv 


and 


Therefore 


From  this. 


- V/) 


> 2.7  X lO'O  X 227  x 1,25  (cm ^/g)  ^ lOOcm^/s^ 

600(2.7  X lO^cmVg) 


a=—  =2xl0  x 1.25  = 25. 


25 


/3  4.7  X 10* 


5.3  X 10  ® cm  =»  53nm . 


Ar(7rP)''^C^>'^. 


227  (tt  X 1,43  X 10  ~^)''^4.18  x 10^  x 2.72  x 10^ 
5.3  X 10  ~*  X 2.7  X lO'®  X 1.25 
-9.7  X 10*s  ^ 


The  time  required  to  form  a bubble  of  radius  /?  — 7 x 10  ^cm  (/  = 


7 X 10 


-5 


5.3  X 10  -* 


= 13.2)  is 


then 


13.2  + In-^  -9.7X  10*/ 


or 

I -4.38X  10  "'^s. 

This  is  the  same  order  of  magnitude  as  that  predicted  by  Seitz  (=10  ” s). 

Clearly  such  small  times  must  be  treated  with  caution.  The  process  is  essentially  an  ex- 
plosion, and  the  difficulty  of  treating  an  untractable  integral-differential  equation  places  a great 
strain  upon  the  analyst  to  quantify  the  growth  of  a bubble  in  such  short  times.  Although  the 
work  of  Forster  and  Zuber  allows  one  to  construct  a graph  of  /(/)  vs  t and  from  it  /(/)  and 
Ht)  needed  for  estimating  the  radiated  sound,  the  limited  conditions  of  validity,  negligible 
heat  input,  simplified  heat  conduction,  mean  value  integration,  etc.)  make  it  unsatisfactory  for 
small  radius  bubbles  and  short  times.  We  adopt  next  a different  approach. 


Approximate  Solution  of  Bubble  Growth  and  Collapse  | 

Lsing  Thermodynamic  Charts  ! 

The  formation  and  growth  of  a bubble  will  be  traced  on  a temperature  vs  entropy  (TS)  i 

chart  for  water  (Fig.  17).  We  imagine  the  process  to  occur  as  follows:  water  at  77°F  (room  tern-  j 

perature),  point  A,  is  heated  nearly  at  constant  volume  (as  if  in  a tiny  steel  tank)  to  620°F,  | 

point  B;  the  heated  spike  then  furnishes  heat  to  an  embryo  along  a constant  temperature  line  | 

BC,  followed  by  a final  expansion  to  ambient  pressure  along  an  unknown  route  CD.  With  ; 

such  a picture  in  mind  we  perform  the  following  calculations.  j 
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I'ig.  17  - Temperature  entropy  chart  for  water 


Embryo  Bubble  Radius  (Critical  Radius) 

At  the  instant  of  explosion  (point  B on  the  TS  chart)  there  is  present  an  embryo  (or  ini- 
tial bubble)  in  the  heated  spike.  We  wish  to  determine  its  size.  The  critical  bubble  size  is  ob- 
tained by  momentarily  balancing  all  pressures  at  the  bubble  wall: 

. I 2<r 

Pv  + Pg  - Poo  + — = 0. 

"cril 

neglecting  the  vapor  pressure  due  to  permanent  gases  ),  we  have 

R = 2(r 

""  P,(T)  -p^- 

Mere  tr  is  the  surface  tension;  estimated  from  Ref.  (21]  to  be  10  dyne/cm  at  T - 600°K.  Tak- 
ing Poo  to  be  15  psi,  we  calculate 

n _ 2x10  , , . „ -7 


(1786  - 15)  6.895  x 10 


— = 1 .6  X 10  cm  — I 6 nil' 
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Mass  and  Volume  of  Water  Heated 


The  mass  of  water  healed  is  determined  by  the  change  in  enthalpy  in  going  from  IV?  to 
620°F.  From  steam  tables  we  find  this  change  to  be  646.7  — 45  = 600  Btu/lb.  For  reasons 
noted  earlier,  we  choose  the  energy  of  the  cosmic  particle  to  be  4 Mev.  Thus  the  mass  of  wa- 
ter heated  (thermal  spike)  is 


m 


= = 4 X lO^eV  X 1.6  X 10  j/eV  ^ _,6 


Ih  600  (Btu/lb)  X 2.32  (J/g)/ (Btu/lb) 


g. 


The  volume  of  heated  water  is  found  from  tiic  specific  volume  of  water  at  620°F,  which’is 
0.0247  ftVib  (X  62.428  = 1.54  cm  Vg).  Thus 

.3 


V^.  =4.6  X 10  ~'^g  X 1.54-^^  =7.1  x 10 


-'‘>cm^ 


g 


whose  radius  is 


= 


7.1  X 10 


-16 


47r 


1/3 


= 5.5  X 10  ^ cm  = 55  nm. 


In  accord  with  the  growth  process  pictured  by  Sette  and  Wanderlingh  this  volume  of  water  is 
considered  to  be  a heat  source  which  contains  a supply  of  heat  (4  MeV)  to  expand  the  embryo 
(process  BCD  on  the  TS  chart.  Fig.  17). 


Final  Bubble  Radius 

From  the  heat  balance  equation  we  have 


AT  = 


niL 


In  infinitesimal  form. 


dr(t)  = 


1 d£(R) 


AT  dE(R) 


dR. 


niL  mL  dR  E dR 

The  minus  sign  means  that  an  increase  in  radius  of  the  bubble  corresponds  to  a decrease  in 
temperature.  As  before 


E(R)  =4TTR^,r*  + nt^  \H'  [TiR)]. 


To  estimate  the  mass  of  vapor,  we  note  that  4 MeV  will  bring  4.6  x 10  g water  to  600°K. 
The  same  4 MeV  will  expand  an  embryo  to  a radius  of  unknown  value  at  point  A on  the  TS 
chart,  where  the  enthalpy  is  1 150  Btu/lb.  Thus  the  mass  of  vapor  is  estimated  to  be 

600 


1150  -45 


X 4.6  X 10 


-16  _ 


= 2.5  X 10 


-16, 


g- 


From  Fig.  18  the  value  of  AW'  of  gassed  water  at  375  K is  approximately  —6.5  kJ/mole. 
Hence  the  energy  required  is 

>7] 


E{R)  =47r(7  X 10  “5)2  140  _ 5.5  x 10^ 

= 5 MeV. 


10 


18 


(2.5  X 10  ■'‘’)  erg 
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where  >r*  = 140  in  from  Ref.  (21).  Here  the  final  temperature  is  taken  to  be  215°F.  Since  this 
energy  is  nearly  the  4 MeV  we  have  postulated,  we  will  accept  the  4 MeV  as  consistent  with  a 
mass  of  2.5  x 10  g vapor. 

Let  us  specify  the  terminal  thermodynamic  point  more  carefully.  The  specific  volume  is 

Y rr  (7  X 10  3 

f",,  ; ^5.74  X 10^-^^  =92ft^/lb. 

^ 2.5  X 10  g 

Thus  point  D on  the  TS  chart  has  h = 1155  Btu/lb,  s = 1 .896,  7'  = 215°F  (—375  K),  and  — 

4,3  psi. 

If  the  water  is  not  gassed,  the  latent  heat  \H''  and  critical  radius  change.  From  a bubble 

at  =7  X 10  cm  is  at  temperature  of  about  410  K (137°C  or  — 280°F),  and  from  Fig. 

18  (Ref.  17)  A//''  is  —4  x 10^  J/mole.  The  total  energy  required  is  still  about  4 MeV,  and  the  ; j 

mass  of  vapor  is  still  2.5  x 10  g.  The  specific  volume  is  also  the  same,  the  enthalpy  is  1185  ] 

Btu/lb,  the  entropy  S — 1.93,  and  the  vapor  pressure  is  4.8  psi.  \ 

.1 

Acoustic  Radiation  Based  on  Rayleigh’s  Formulas  I 

Since  the  derivation  of  a history  of  wall  expansion  (namely  Rit))  is  fraught  with  \ 

difficulties,  we  will  resort  to  a modified  Rayleigh  theory  to  calculate  the  acoustic  effect  of  grow- 
ing and  collapsing  bubbles.  A convenient  summary  of  this  modified  theory  is  found  in  Ref. 

(221. 

We  first  simplify  the  picture  and  assume  there  is  an  average  driving  pressure  throughout 
the  growth  period  which  performs  work  in  expanding  the  bubble  from  R — O \o  R = 1 

X 10  cm.  We  neglect  compression  of  the  medium.  Thus,  for  a volume  change  A I , and  ! 

100"/)  heat-to-work  efficiency, 
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^</v  ^ ^ ^ inpul 


or 


= 


4 X 10^  X 1.6  X 10  -1^  era 
j 7r(70  X 10 


= 4.45  X 10*"  dyne/cm 


Thus  although  the  initial  pressure  is  1786  psi  abs  and  the  final  pressure  is  about  4 psi,  the 
average  pressure  is  about  4-1/2  atmospheres  (say  67  psi).  The  potential  energy  of  the  fully  ex- 
panded bubble  is  then 


P.E.  = « A f =6.4  X 10 


erg. 


According  to  Rayleigh’s  theory  the  wall  velocity  is  approximately  constant  during  most  of  the 
growth  phase.  Thus  [23] 


4.45  X lO^dyne/cm^  cm'* 
1 dyne  ■ s^ 


= 17  m/s. 


(Note  that  the  density  of  the  water  is  taken  at  room  temperature.)  Since  this  wall  velocity  is 
well  below  the  velocity  of  sound  (-1.5  km/s),  compressibility  effects  can  be  neglected  The 
average  time  required  to  expand  the  bubble  to  its  final  radius  is 


0.7  X 10 


m 


0.04  /.IS. 


17  m/s 

The  average  acoL'<-tic  power  radiated  over  the  period  of  the  growth  phase  is  obtained  by 
legrating  over  the  square  of  the  acceleration: 


=■ 


4tt  C, 


f Pit)  dt.  ii  =i,-!t(2RR^  + R-ii 


/ 0 


). 


For  constant  velocity  this  reduces  to  [24] 

E,... 


P.E. 


= 1 

^ y 3 y 


4.45  X 10^ 

1 (1.5  X 10^)2 


0.03. 


Thus  only  about  3%  of  the  input  energy  is  radiated  as  sound.  For  an  input  of  4 MeV  the 
acoustic  power  radiated  in  0.04  ix%  is 


=0.03  X 4 MeV  x 0.16  pJ/MeV  = 0.02  pJ  = 2 x 10 


erg. 


The  average  acoustic  pressure  at  1 m from  the  center  of  the  bubble  is 

1/2  ^ ..  / ^a<  ^ / 

\ AM7r/?2  y 


< p^> 


2 X 10  2 X 1 X 1.5  X 10^ 


4 X 10  4rr  (10^)2 


or 


Pav  “ dynes/cm  2 at  1 m,  transient,  from  one  bubble  due  to  4-MeV  cosmic-particle  interac- 
tion product. 
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Resonant  Frequency,  Steady-State  Radiation 


The  pressure  calculated  above  is  the  average  (shock  wave)  pressure  transient  over  the 
period  of  bubble  formation.  Subsequent  to  formation  the  bubble  will  execute  quasi-harmonic 
vibrations  at  the  resonant  frequency  .//i;,  where  [25] 


Jr 


_J 

2rr/?„ 


in  which  is  the  nominal  bubble  radius,  y = C^IC\.  of  the  gaseous  contents  of  the  bubble, 
Pp,.,  is  the  external  pressure  on  the  bubble,  and  p,  is  the  density  of  the  water  "at  infinity  ."  For 
air  = 0.24,  Q = 0.17,  and  y = 1.41;  for  water  = 0.35,  = 0.27,  = 1 3 at 

room  temperature.  For  ease  in  computation  we  take  y = 4/3  for  both  gases  in  the  cavity.  The 
external  pressure  is 


P = P 

* C.\l  ' a 


2<r  ( T) 

R.. 


Let  us  make  this  calculation  on  the  assumption  that  the  final  bubble  temperature  is  410  K 
(137°C  = 280°F);  we  then  estimate  (r  =40  dyne/cm  (21).  Sei’ing  P^  = 1 atmosphere,  we  ob- 
tain for  a 0.7//m  bubble, 

=1.03  X 10*  = 22  X 10*dyne/cm2 

7 X 10  -5 

Thus  the  resonant  frequency  of  the  bubble  is 


Jr 


! 

277  X 7 X 10 


= 7 MH£. 


To  calculate  the  acoustic  radiation  due  to  resonance  of  the  bubble,  we  must  account  for 
energy  losses.  We  assumed  earlier  that  the  potential  energy  of  the  bubble  at  radius 
R = 0.7  /ixm  is  exactly  the  energy  of  deposition  (4  MeV).  Now  let  y represent  the  fraction  of 
energy  input  lost  in  shock-wave  radiation  both  in  the  growth  and  collapse  phase.  Although  we 
have  not  calculated  the  collapse  phase,  we  assume  it  radiates  the  same  as  the  growth  phase. 
Thus  the  energy  available  for  energizing  bubble  vibrations  is  = (1  - v)  = (1  - v )4 

MeV.  We  use  this  energy  to  determine  the  volume  displaced  in  going  from  /?=/?„= 

= 0.7  /am  (say)  to  some  minimum  radius  R„  — a during  vibration.  The  potential  energy  at 
minimum  radius  is  seen  to  be 


so  that 


P.P.mm  = 67ryP,,^,P„«^ 


,/Zj7  “ ■ 

Y (’77yP,,,P„ 


I 


1 

j 


I 


i 
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I 

K- 


1 


i 

► 


i 

I 


where  is  the  potential  energy  ar  radius  R.  From  the  previous  calculation  we  estimate  >■  to 
be  0.06,  so  that  =0.94  x 4 >:  10**  x 1.6  x 10  =6  x 10  erg.  Hence 


6 X 10 


= 0.39  Mm. 


2.2  X 10*’  X 7 X 10 


ir  we  assume  sinusoidal  radial  motion  during  vibration,  so  that  at  the  resonant  frequency 


then 


R = R„  + (1  sin 
R = In/ffU  cos 


The  volume  velocity  is  then 

S = {AttR^,  )2itI i^a  cos 

Hence  the  steady-state  amplitude  of  radiated  (acoustic)  pressure  will  be  (at  1 m) 


or 


l/’,J 


0,1,  f 

4nR 


(<...S) 


(27r  X 6.7  X ) (7  x 10  (3.9  x 10  ) 

10^ 


I =3.4  dynes/cm  rms  steady  state  (maximum) 

at  1 meter  at  resonant  frequency  //;  = 6.7  MHz 
from  1 bubble  excited  by  0.94  x 4 MeV. 


Since  the  radiation  has  a finite  band  of  frequencies  because  the  radial  motion  is  not  exactly 
periodic,  we  will  calculate  the  acoustic  effect  at  25  kHz.  That  is,  roughly. 


1 n 1 

= 3.4  X 

o 

X 

2 

= 4.7  X 

10  'dynes/cm^ 

1 1 

6.7  X lO^* 

rms  steady  state  at  25  kHz  at  1 
excited  by  0.94  x 4 Me\'. 

meter  from  1 bubble 

This  result  is  based  on  acoustic  pressure  being  proportional  to  radial  acceleration  of  the  bubble 
wall,  hence  proportional  to  (frequency  ).^ 
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Extrapolation  to  Higher  Energies 

Let  us  assume  that  a single  cosmic  particle  enters  the  ocean  with  an  energy  ol  10  eV.  If 
every  4-MeV  interaction  collision  generated  one  bubble,  we  could  expect  2.5  x 10^  bubbles  to 
form  (at  the  maximum).  Further,  if  all  bubbles  radiated  coherently,  then  the  sieady-siate 
acoustic  pressure  measured  at  25  kHz  would  be 


= 4.7  X 10  “5  X 2.5  X 10’  = 1.2  X 10  ’ dynes/cm  ’ 
meter  at  25  kHz  due  to  a 10  '‘’-eV  particle  (steady  stale). 


At  100  m this  pressure  reduces  to 

\p„^\  = 12  dynes/cm 
This  is  above  ambient  noise  (at  25  kHz). 

Collapse  Pressure  and  Acoustic  Radiation 


The  initial  calculation  made  above  was  based  on  an  isothermal  expansion  BC  to  a final 
state  flat  280°F  and  about  4 psi,  the  last  part  of  the  expansion  being  unknown.  To  estimate 
collapse  pressures,  we  require  a more  expicit  statement  of  the  final  stage  of  expansion.  From 
Ref.  17  it  can  be  inferred  (though  it  is  not  so  staled  by  them)  that  the  final  expansion  is  adia- 
batic. Hence  we  hypothesize  the  path  on  the  TSchart  (Fig.  17)  to  be  ABCDBf. 


Let  us  calculate  the  final  radius  at  point  E.  We  calculated  previously  that  4.6  x 10  g 
of  water  are  heated  on  the  thermal  spike  to  620°F,  requiring  an  input  of  4 MeV.  We  now  allow 
this  heated  spike  to  transfer  heat  to  the  embryo  bubble,  expanding  it  isothermally  from  flto  ( 
to  £,  the  point  £ being  at  620°F  and  atmospheric  pressure,  with  a specific  volume  of  43  6 ft  ’/lb 
(x  6.24  = 2.72  X 10  ’ cm ’/g),  enthalpy  of  1344  Btu/lb,  and  entropy  of  1.982.  We  next  ima- 
gine that  the  final  expansion  from  £ to  £ takes  place  at  constant  entropy  to  a temperature  of 
280°F.  At  this  point  the  pressure  Py  is  2.1  psi,  the  specific  volume  is  192  ft  ’/lb,  and  the  enthal- 
py is  1159  Btu/lb.  To  bring  water  to  superheat  at  £ requires  an  addition  of  1344.5  — 45 
1300  Btu/lb.  We  estimate  the  number  of  grams  of  superheated  vapor  to  be 

"'v  = -^r  X 4.6  X 10  =2.12  x 10  g. 

1 300 

Since  the  specific  volume  is  192  x 62.4  « 1.198  x lO'^cm’/g.  the  volume  of  superheated 
vapor  is 
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This  is  close,  but  not  quite  the  assumed  final  value  of  0.7  moi.  To  obtain  a more  consistent  final 
point  we  require  that  the  adiabatic  expansion  arrive  at  a radius  of  0.7  /im,  hence  at  a specific 
volume  of 


y 77(7  X 10-5)^ 

2.12  X 10 


cm  ^ 
g 


= 108  ft^/lb. 


= 6.8  X 10^  cm  ^/g 


The  terminal  point  (from  steam  tables)  is  then  350°F,  enthalpy  is  1209  Btu/lb,  entropy  is 
1.982,  and  pressure  is  4.5  psi. 


After  expansion  the  forces  acting  on  the  maximum  bubble  are  unbalanced.  The  net  col- 
lapse pressure  will  now  be  calculated  This  is 

n I 2(7 

/>,  - P^. 

at  T = 350°F  (=175°C)  we  estimate  it  = 35  dyne/cm,  so  that 

P,r  = = 10  dyne/cm 

7 X 10-5 

Also,  at  a temperature  of  350°F  ( = 458  K)  we  estimate  the  dissolved  gases  to  have  a pressure 
of  2.7  X 10^  dynes/cm  ^ (see  below  for  this  estimate).  Since  = 4.5  psi  ( x 6.895  x lO"*  = 
3.1  X 10  5dynes/cm  ^),  and  since  = 1 atmosphere, 

P,,,,  = 1 X lO**  -f  1 X lO^"  -3.1  X 105  -2.7  X 105 

= 1.69  X 10*dynes/cm 

Let  us  asume  that  during  collapse  the  vapor  is  not  instantly  reabsorbed,  so  that  it  can  act  as  a 
buffer  to  total  disappearance  of  the  cavity.  Then  the  ratio  of  the  net  collapse  pressure  to  the  va- 
por pressure  (called  P/Q)  is 


Q 


1.69  X 10^ 
3.1  X 105 


= 5.44,  or 


SI 

p 


= 0.18. 


These  ratios  will  enable  us  to  calculate  bubble  dynamics  on  collapse.  The  minimum  radius  of 
collapse  is  obtained  by  setting  (261 


R = /? 
^niin  ^max 


Q 1 ' 5(>  -1  ) 

II  \ 

Piy  - 1 ) 1 

^ /•(>  - 1 ) I 

-I /3(y  -I  ) 


This  formula  is  due  to  Neppiras  and  Noltingk  (27). 


Assuming  y =4/3,  we  obtain 


0.18  X 3 
1 -I-  0.18  X 3 


= 0.35 


S«ty 


This  paramciei,  plus  knowledge  of  P/Q,  allows  us  to  calculate  the  peak  radiated  shock  wave  on 
collapse.  From  Ref.  [28]  at  1 m 


P R 

+ _ ^max 

D 

^max 

1 

1 -4 

^min 

3 

"max  -j 

3 r 

^min 

D 

^max 

1.69  X 10*  dyne/cm  ^ 

' 10  -5cm 

(3)^ 

1 

3 

X — 

10^ 

1 
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or 


^niax  — ^ dyne/cm  ^ shock  amplitude  at  1 m during  collapse  of  a 

0.7  m bubble  based  on  assuming  the  vapor  acts  as  a buffer. 


The  shock  power  radiated  is  (roughly)  given  by  [29] 


Ea.  1 

/ 2F,„„ 

D 

''max 

P.E.  3 I 

/ h>fC} 

^niin 

■i/ 

= 0.012. 


2 X 1.69  X lO^* 

3 X 1 X (1.5  X 10^)2 


We  estimated  earlier  that  the  acoustic  output  on  collapse  was  about  the  same  as  upon  growth 
(3‘Ki  each).  Mere  the  figure  on  collape  is  1%.  Of  course  some  of  the  vapor  will  be  reabsorbed  on 
collapse.  Hence  I’Ki  is  too  pessimistic.  If  all  the  vapor  returned  to  the  liquid  (water)  upon  col- 
lapse, and  left  only  permanent  gas  pressure,  then  it  is  simple  to  calculate  that  the  peak  shock 
wave  pressure  will  be  = 1.5  x lO^dynes/cm^  at  1 m from  one  bubble.  This  is  overly  op- 

timistic. 


Calculation  of  Gas  Pressure 


The  partial  pressure  of  the  permanent  gases  in  the  bubble  is  an  essential  parameter  in 
calculating  the  shock-wave  magnitude.  However  it  is  difficult  to  calculate  it.  The  gas  content  of 
the  ocean  is  reported  in  books  on  oceanography  in  the  following  way.  First,  whatever  the  loca- 
tion of  the  water  volume,  it  is  assumed  that  at  some  time  it  had  been  at  the  surface  of  the  oce- 
an and  in  equilibrium  with  the  air  [30].  Thus  regardless  of  depth  the  water  is  assumed  nearly 
saturated  (at  NTP  = normal  temperature  and  pressure).  Hence  gas  content  is  reported  in  sa- 
turation values  (say  nil/l)  in  equilibrium  with  a normal  dry  atmosphere.  From  page  188  of  Ref. 
[21]  we  see  that  at  0°C  and  35%  salinity  the  saturation  volume  of  oxygen  in  seawater  is  about  8 
/»/// and  of  nitrogen  is  about  14  »;///.  We  next  suppose  this  dissolved  gas  is  in  the  microbubble 
de.scribed  above,  and  we  desire  to  calculate  the  partial  pressure  of  the  dissolved  gases  at  the  in- 
terface. According  to  Henry's  law  the  concentration  m of  a gas  in  a liquid  is  related  to  the  par- 
tial pressure  of  the  gas  as  m = C\p,  where  C\  is  the  coefficient  of  saturation.  If  m is  expressed 
in  units  of  milliliters  per  liter  and  p is  expressed  in  torr  (760  torr  = 1 atmos  = 1 
X 10^  dynes/cm  ^),  then  at  0°C  (273  K)  and  35%  salinity  the  value  of  C\  is  38  for  oxygen 
(Oj)  and  14  for  nitrogen  (Nj)  [Table  41,  p.  191,  Ref.  [30]].  Thus  the  partial  pressures  are 

■T60~'  “=2.8  X 10  ^ dynes/cm 


and 


I 


11  X 1 X 10^ 

14  760 


= 1 .3  X 10^  dynes/cm 


The  total  gas  pressure  is  therefore  about  1.6  x 10^  at  0°C.  However  the  temperature  of  the 
bubble  at  maximum  radius  is  taken  to  be  458  K.  Hence  the  estimated  gas  pressure  is  458/273 
X 1,6  X 10^=2  6 X 10^  dynes/cm  ^ 
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Effect  of  Depth  of  Particle  in  the  Ocean 


Let  us  again  consider  a 4-MeV  particle  heating  a mass  of  4.6  x 10  g water  in  a ther- 
mal spike  to  620°F,  at  which  temperature  the  specific  volume  is  0.0247  fi^/lb  (=  1.54  cm  Vg). 
As  before  this  water  acts  as  a heat  source  and  feeds  an  embryo  that  expands  (say)  to  ambient 
pressure.  Let  ambient  be  1600  psi  (say  about  3500  ft,  or  1 km,  of  water).  At  the  end  of  ex- 
pansion the  pressure  is  so  great  that  the  specific  volume  is  only  0.272  ft^/lb  (x  62.4  = 16.97 
cm^/g)  the  enthalpy  is  1187  Btu/lb,  and  the  entropy  is  1.348.  Since  only  600  Btu/lb  is  avail- 
able, the  mass  of  the  heated  vapor  is 

600 


/»/„  = 


1187  -45 

Thus  the  volume  of  superheated  vapor  is 


X 4.6  X 10 


-16  _ 


= 2.4  X 10 


-16 


g- 


y.  = 2.4  X 10  "'^g  X 16.97 


cm  ■ 


= 4.1  X 10 


cm 


whose  radius  is 


= 


4.1  X 10  “'5  X 


I 

47r 


1/3 


= 9.9  X 10  * cm  ( say  0.1  /urn  ). 


We  compare  this  expansion  to  1600  psi  with  the  earlier  calculated  expansion  to  15  psi  (at 
620°F).  At  the  latter  point  the  specific  volume  is  2.72  x 10^  cm  ^/g,  from  which  the  volume 
of  heated  vapor  is 


F,  =2.4  X 10  "'^g  X 2.72  x 10^-^^  =6.5  x 10~'^cm^ 

g 


whose  radius  is 


'15/.V 


6.5  X 10 


-13 


47r 


1/3 


= 5.4  X 10  ^ cm  =0.54  urn. 


Thus  the  rcsultani  bubble  at  1600  psi  is  about  5-1/2  times  smaller  than  at  15  psi  for  the  same 
energy.  However  it  is  larger  than  critical  (=  1 nm);  hence  it  can  grow  by  further  adiabatic  ex- 
pansion. But  Settc  and  Wanderlingh  show  [Fig.  7]  that  the  final  expansion  (after  isothermal 
expansion)  is  small.  We  have  calculated  that  in  the  case  of  adiabatic  expansion  from  620‘’F  at 
15  psi  to  350°F  at  4.5  psi  the  change  in  radis  is  about  2/7  or  30"/:.  Thus  we  estimate  the  final 
radius  at  the  depth  of  1 km  to  be  about  0.13  still  a factor  of  about  7 smaller  than  a bubble 
near  the  ocean  surface.  This  smaller  radius  bubble  will  reduce  the  acoustic  eflfect.  For  exam- 
ple the  peak  shock  wave  on  collapse  will  be  about  7 times  smaller,  assuming  the  ratio  of 
^min' is  about  the  same  as  in  the  earlier  case.  Greater  depths  will  even  have  a more 
drastic  effect,  primarily  due  to  reduction  of  the  maximum  bubble  that  can  form  with  the  given 
input  energy.  Of  course,  if  higher  energy  particles  participate  in  bubble  formation,  a larger 
acoustic  output  can  be  expected. 


JET  MODELS  AND  OTHER  MODELS 
Jet  Models 

After  a cosmic  particle  creates  a shower  of  interaction  products  in  a liquid,  the  shower  it- 
self can  be  imagined  to  be  the  continuation  of  a jet  of  equal  diameter  entering  the  liquid  from 
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the  outside.  This  jet  creates  noise  as  it  traverses  the  liquid.  Several  mechanisms  may  be  con- 
sidered. The  principal  mechanism  is  injection  of  mass  of  jet,  hence  monopole  radiation  is  to  be 
expected.  However,  if  the  flow  of  the  jet  meets  obstacles  (fluid  molecules),  the  expected  radia- 
tion is  dipole.  If  the  jet  creates  thermal  (or  velocity)  turbulence,  the  radiation  will  be  quadru- 
pole.  Our  objective  will  be  to  calculate  the  noise  power  due  to  these  three  mechanisms. 

By  dimensional  analysis  one  can  establish  relations  between  the  parameters  of  jet  noise. 
These  are.  the  density  of  the  jet  material,  the  speed  V of  the  jet  flow,  the  pertinent  length  L 
and  the  mach  number  M = V/c\  where  c is  the  speed  of  sound  in  the  liquid.  For  jets  of  finite 
cross-sectional  area  S one  chooses  L such  that  = S.  The  relations  between  the  noise  power 
W' generated  and  these  parameters  are  as  follows: 

monopole  radiation;  (F,,,  oc  M 

dipole  radiation;  a.  y^ 


quadrupole  radiation: 


a p,  /,2  yi  a/5 


The  constants  of  proportionality  depend  on  the  process.  Various  estimates  of  their  magnitude 
have  been  calculated.  For  example,  if  the  mechanical  stream  of  the  jet  has  a kinetic  power 

^imxh  (units  of  w^:  kg/s). 

then  the  acoustic  power  generated  is  roughly 

^ IF ^ = ,0. 


A/55 

!Ll 

rJi  ,420- 

ft 

a/5-5 

10.830’ 

ft 

A/55 

JIl 

120,470’ 

ft 

5F  = 

M 

provided  < ^^,,^,^11  (Ref.  31).  Since  = p,  SF  =(1^  F,  we  see  that 

M",/  = const,  M 
= const.  M-y, . 

The  constant  in  each  case  matches  the  constant  in  the  expression  for  IF^,. 

The  calculation  of  the  effective  mach  number  M of  the  fluid  jet  requires  considerable 
care.  Let  the  imagined  cosmic-particle  jet  have  a unit  flux  of  energy  (per  cross-sectional  area) 
of  magnitude  H',.  Upon  interaction  with  the  liquid  the  particle  jet  is  converted  to  a liquid  jet 
whose  power  is  IF^.  The  interaction  itself  depends  on  the  capture  area  presented  by  the  fluid 
to  the  incident  particle  stream.  The  ratio  of  powers  IF,  to  IF,,  is  designated  y,  where 

m,  y}n 

^ ~ nt,  F,/2  ’ 

in  which  and  w,  represent  ma.ss  flow.  Assume  next  that  the  mass  flows  are  proportional  to 
mass  density.  Then 


or  F,  = 
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This  is  the  initial  velocity  of  the  fluid.  Thus  the  effective  mach  number  of  the  jet  is. 

Noting  that  ff',,,,,,/,  = y we  see  that  the  acoustic  power  generated  is 


where 


= 1420 

>>1 

= 10 

^'l 

= 10,830 

Rj 

f’l 

fh 

= 1 

= 120,470 

yi 

= 0.1 

The  coefficient  of  momentum  transfer  y is  roughly  the  ratio  of  the  scattering  cross  section  of 
the  incident  particle  to  that  of  the  fluid  particle.  The  noise  power  generated  is  seen  to  be  pro- 
portional to  y 


Transient  Radiation  From  Sources  in  Motion 

When  a cosmic  particle  is  decelerated  by  a liquid  the  liquid  itself  is  accelerated.  Such  an 
acceleration  constitutes  an  acoustic  source,  and  an  acoustic  pressure  (shock  wave)  is  radiated 
outward.  A simple  theory  constructed  on  this  picture  is  as  follows. 

Let  the  volume  flux  of  liquid  set  into  motion  be  S(/»(Vs),  and  let  the  corre'^ponding  mass 
flux  be  ii/S  = </(/)  (units:  N • s/m).  The  motion  of  the  cosmic  particle  is  equavalent  to  the 
motion  of  </(r).  The  latter  constitues  source  distribution  density  0(r,r)  (units:  N s/m*)  of 
form, 

Q(r.r)  - yf)  S(r)S(x). 

which  radiates  shock  waves  according  to  the  formula 

^2/,  _ J.  li£  = _A^(,)S(r  - f'r)  fi{,v)  fi(.v). 

The  solution  of  this  equation  is  obtained  by  a coordinate  transformation  (4,  pp.  721  ffi.  When 
the  velocity  V of  the  equivalent  source  is  supersonic  (as  we  assume),  the  solution  can  be  ex- 
pressed as  function  of  distances  R R and  angles  « explained  below;  thus 


■f  = --^^(,)S(-  - t',)  fi(.v). 

nt 


p(R.H.l)  -f  

Att  R (A/cos  W — 1 ) ^ An  R (m  cos  H — 1 ) ^ 

_ </(>  ~ /?  ^/c)  (A/  - cos  0 ) y ^ q{i  - R ~/c)  (M  - cos  R ~)  f 

An  (R*-)^  (A/ cos  « + - 1 An  (R  -)^  (A/ cos  H “ - I )^  ’ 
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I ig.  19  — Transient  radiation  from  a source  in  motion 


in  which  q'  is  the  derivative  with  respect  to  the  argument.  Distances  and  angles  are  given  by 
Fig.  19.  In  this  figure  the  source  S is  moving  in  the  direction  r with  supersonic  velocity  (' 
(mach  number  M > \).  It  is  the  origin  of  a cone  of  shock  waves  (A,  B).  Any  point  P inside 
the  cone  will  hear  the  pressure  p given  above  at  time  /.  This  pressure  will  consist  of  radiation 
from  two  fictitious  origins,  E and  £_,  with  at  a distance  R from  P,  and  E _ at  dis- 
tance R _ from  P (E  _ is  in  front  of  /’and  £4.  is  behind).  The  distances  are  given  by 


The  half  angle  of  the  cone  is  w = sin  (W  ).  Outside  the  cone  at  time  1 there  is  silence. 
The  directionality  inherent  in  the  formula  for  pressure  p must  be  obtained  by  numerical 
evaluation.  If  the  time  dependence  is  harmonic,  ftiat  is,  q'(t  — R/c)  = cos 

(o„(r  — R/c)  then  the  fields  due  to  £4.  and  £_  interfere,  so  that  the  reserved  sound  pressure 
oscillates  with  time  (beating  phenomenon).  The  amplitude  spectrum  of  the  radiated  pressure 
can  be  obtained  by  Fourier  transformation: 


p(RM,(o ) 


' p(r,P,l)dl. 


Another  model  of  cosmic  particles  as  travelling  sources  may  have  relevance  to  acoustic- 
pulse  formation.  In  this  picture  the  interaction  products  bounce  back  and  forth  due  to  succes- 
sive collisions.  Thus  the  advance  of  the  cosmic  particle  is  imagined  to  be  that  of  a dipole 
source  travelling  with  supersonic  speed.  A formulation  of  the  pressure  field  developed  is 
presented  for  convenience. 


The  bouncing  byproducts  exert  a force  £,  (per  unit  volume)  conceived  to  be 


£,  =./,  (r)  fi(r  - Vt)  fi(.v)  fi(>  ), 

in  which  ,/■  (r)  is  an  impulse  force  (units:  N).  To  apply  this  force  to  the  pressure  equation,  one 
lets  p =div  A,  and  solves  for  A;  the  governing  equation  is 
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1 d ^ 4 

V^/1,  =./;(')  8 (-'  - P'r)  8(jf)  8(>). 

b I 

whose  solution  is 

^ ^ _ f,  - fi  ^/c)  J]  0 - R -/c) 

477  /?,+  477  /?,- 

Ultimately  the  field  pressure  is  [4  page  721fn 

f/U  - R ■'■/c  icos  « f- {i  - R +/[  ) (M  - cos  I)  +) 

477  f (A/ cos  « —1)^  477  (/?■'')  ^ (A/ cos  « —1)^ 

J]  'U  - R 7c)  cos  a - lAi  - R 7c)  (A/  - cos  H ) 

+ —L ^ 

477  /?  c (A/  COS  « - 1 ) 4ir  (R  ) ^ (A/  cos  « ~ - 1 ) ^ 

This  dipole  case  is  conceptually  related  to  the  monopole  travelling  source;  if  one  sets  q V =./, 
so  that  the  oscillating  force  is  numerically  equal  to  the  monopole  mass  flux  traveling  at  velocity 
T,  it  is  then  possible  to  interchange  the  formulas.  In  particular  one  sets  q'  and  then  the 

derivations  of  the  formulas  are  identical. 

CONCLUSION 

The  generation  of  sound  by  sources  is  well  understood  when  the  sources  can  be  categor- 
ized as  monopole,  dipole,  quadrupole,  etc.  In  the  particular  case  of  the  noise  pulses  generated 
by  cosmic  particles  in  the  ocean  the  greatest  emphasis  in  modeling  to  date  is  on  monopole 
sources  related  to  the  time-varying  application  of  heat.  Two  varieties  of  heat  models  have  been 
surveyed  and  developed  in  this  report:  the  thermoelastic  model,  in  which  a deposition  of  heat 
into  the  liquid  results  in  an  elastic  expansion  which  radiates  outward  as  a shock  wave,  and  is 
the  microbubble  model,  in  which  the  deposition  of  heat  "boils"  the  liquid  locally,  creating  an 
expanding  (and  collapsing)  bubble,  which  radiates  sound  as  shock  waves.  Other  models  are 
possible,  such  as  a jet  noise  model,  a moving-source  model,  and  an  explosion  model.  The  first 
two  have  been  noted,  but  no  detailed  calculations  given.  An  explosion  model  is  discussed  in 
Appendix  A.  Numerical  calculations  of  several  heat  models  have  been  made  in  the  text  and 
are  further  amplified  in  Appendixes  B and  C.  As  expected,  they  yield  different  (even  radically 
different)  answers  for  the  same  input  parameters.  These  differences  are  attributed  to  different 
conceptions  of  what  the  precise  physical  event  is  in  the  liquid.  It  is  concluded  that  further  pro- 
gress in  modelling  must  await  renewed  effort  in  experiment  work. 
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Appendix  A 

RADIATION  MODEL  BASED  ON  THE  THEORY  OF  EXPLOSIVE  SOURCES 

In  the  theory  of  explosive  acoustic  sources  the  potential  energy  of  the  explosion  bubble 
(iFg)  is  only  a small  percentage  of  the  intrinsic  energy  of  the  source  that  is. 

iCg  = A"  iCj,  where  K is  0.10  or  less.  In  the  well-known  theory  of  Willis  the  initial  radius  of  a 
bubble  due  to  potential  energy  IFg  at  hydrostatic  pressure  Pis 

I 1/3 

■ rr  n I ■ ^ A IV  I 

R.  = 


1/3 

3 A IF  , 

AttP 

AnP 

and  the  first  period  of  the  bubble  oscillation  from  the  Willis-Rayleigh  theory  is 

. j 1/2 

Tfl  = 2 X 0.914  ■ ■ 


1/3 

^>l 

AttP 

P 

Two  types  of  radiation  propagate  from  the  bubble;  a shock  wave,  and  tnen  a wave  due  lo  bub- 
ble oscillation.  The  peak  pressure  developed  by  the  shock  wave  has  been  empirically  found  to 
be 


/’max 


W 


1/3 


R 


m which  IF  is  the  weight  of  explosive  (its  internal  energy),  R .s  the  distance  of  measurement, 
and  k,  a are  empirical  constants.  Many  experiments  show  that  « is  approximately  unity.  Thus 
/max  * ^ [Ala].  To  include  the  effect  of  the  explosion  pressure  subsequent  to  the  passage 
of  the  peak  pressure,  one  uses  the  time  integral  of  the  pressure,  called  the  impulse,  defined  as 

I 

/(/)  = r Pii)  dr. 

*0 

Empirically  the  impulse  is  given  by  [Alb] 

/ = / 

where  / and  /3  are  constants 


IF 


1/3 


A third  characterization  of  the  shock  wave  is  the  energy  flux  E,  per  unit  area  of  a fixed 
surface  normal  to  the  directions  of  propagation.  This  is  given  by  [Ale] 


E,  = m 


W 


1/3 


R 


where  m and  y are  constants 


To  use  these  empirical  formulas,  one  requires  a knowledge  of  the  constants.  They  will 
depend  on  the  type  of  explosive  and  the  units.  For  example,  let  the  explosive  be  TNT  (densi- 
ty 1.52)  and  the  units  be  English:  pressure  in  psi,  distance  in  feet,  weight  in  pounds  (mass), 
and  energy  density  in  inch-pounds  of  force  per  square  inch.  Then  [Aid] 
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i 

i 

[ 


and 


= 2.16  X lO'* 


lb 


/ =1.46 


H.1/3 


R ft 
0.89 


R 


1.13 

(units:  psi ), 
(units:  psi  • s ) 


E = 2.41  X 10-’  W 


1/3 


^^l/3 


R 


2.05 


units:lb  — in./in. ^ 


Here  the  integration  time  is  taken  (somewhat  arbitrarily)  as  6.7W,  where  W is  the  time  constant 
of  the  initial  high-pressure  region. 


To  use  these  formulas  for  the  case  of  high-energy  particles  we  must  relate  the  particle  en- 
ergy to  the  weight  of  TNT  that  is  energetically  equivalent.  Only  approximations  are  possible, 
and  these  depend  on  experience  factors.  Since  the  energy  release  of  TNT  depends  on  depth, 
we  will  take  30  ft,  or  9 m,  to  be  the  depth  in  the  following  calculation. 


We  assume  from  experience  lA2,  p.  17]  that  1 lb  of  TNT  will  release  roughly  1.5  x 10^ 
ft-lb  of  energy.  (Another  convenient  estimate  is:  1 g of  explosive  liberates  1 kilocalorie  = 
4.18  kJ  =2.61  X 10^^  eV,  or  1.19  x 10^^  eV/lb  or  1.4  x 10^  ft -lb.)  Converting  to  more  con- 
venient units,  we  have  for  the  conversion  factor  E^^/ 

^ _ I.S  X 10<>ft.  Ib/lb  X 1.3558  J/ft.  lb 

” 453.6  g/lb  1.6  X 10  J/eV 

= 2.802  X 10^2  j 271  x 10“eV/lb. 

g 

For  example  a high-energy  particle  with  £„  = 10 eV  is  equivalent  to 

1 0 pV 

— If =7.87  X 10  “">16  of  TNT 

1.27  X iu25eV/lb 


or 


10 '^eV 


= .0356  #<g  of  TNT. 


2.802  X 10^2  eV/g 

In  conformity  with  our  calculations  in  the  section  on  microbubbles,  we  apply  the  above  formu- 
las to  the  case  of  a 4-MeV  particle.  The  equivalent  weight  of  TNT  of  such  a particle  is 

4 X 


W = 


10^  eV 
25  , 


1.271  X lO^^eV/lb 


= 3.15  X 10  '’’lb  of  TNT. 


We  indicated  in  the  main  text  that  the  conversion  efficiency  from  particle  power  to  acoustic 
power  is  roughly  3%  on  expansion  of  the  bubble  and  1%  on  contraction.  Let  us  arbitrarily  as- 
sume a 2%  energy  conversion.  Thus  the  actual  energy  available  to  form  a bubble  is  assumed  to 
be 


fF  = 3.15  X 10  X 0.02  = 6.3  x 10  lb  of  TNT. 


Thus,  applying  the  formula  for  the  peak  shock  pressure  (at  1 m = 3.28  ft),  we  obtain 
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J (6  3 X 10  -21)  1/3  ‘ 

=2.16  X lO'*  = 1.388  X 10  psi 

= 9.58  dynes/cm ^ (peak). 

Similarly,  we  can  calculate  the  impulse  at  I m: 


/ = 1.46  X (6.3  X 10  -21)  1/3 


(6.3  X 10  -21)1/3 


— 9.52  X 10  psi -s  or  6.57  x 10  ^ d>  nes/cm  ^ s. 

If  the  dominant  frequency  is  25  kHz,  we  assume  the  integration  time  to  be  1/25,000  = 40ms. 
The  averaf;e  pressure  radiated  is  then 

„ - J - 6.57  X 10  dynes/cm^  s dvne  , 

/’ov  7^  TTZI I 64  10  ^ - at  I meter. 

40  X 10  s cm‘ 

We  note  that  both  and  are  of  the  sam.'  order  o'  magnitude  as  previously  calculated  in 
the  main  text.  However  this  choice  of  integration  time  is  quite  arbitrary.  Cole  [Ale,  p.  239) 
takes  the  integration  time  to  be  r = 6.7«  where  «,  is  the  time  constant  of  exponential  decay  of 
the  main  shock  pulse  (as  if  the  pulse  had  the  form  p = />,„  exp  (-//«)).  A plot  of  reduced  time 
constant  tt/lf'’/^  vs  [Alf,  p.  240]  is  out  of  the  range  of  our  numerical  work.  However 

we  can  try  an  estimate.  We  note  that  «/ir'  ^ rises  0.05  ms  sec  per  falling  decade  of  W^''^IR. 
For  a 4-MeV  particle  whose  effective  TNT  equivalent  is  roughly  6.3  x 10  lb,  the  value  of 
W^'^/R  at  /?  = 1 m is 


(6.3  X 10  -21  lb)  1/3  1 8469  x 10 


= 5.6  X 10 


3.28  ft  3.28 

This  is  roughly  10^  or  7 decades  lower  than  W^  '^/R  = 1.  Hence  by  extrapolation 
= 0.05  X 10  + (0.05  X 10  ~^)  10^  =5  x 10^. 

The  time  constant  is  therefore 

M = 5 X 10^  = 5 X 10^  X (6.31  x 10  = 92.3  ms. 

Our  estimated  integration  time  is 

I = 6.7  X 0 = 6.7  X 92.3  ps  — 0.62  ms. 

The  average  pressure  radiated  in  this  time  is  estimated  to  be 

— f - 6.57  X 10  dynes/cm  ^ s , -s  . , i 

Pm  ~ =1.06  X 10  dynes/cm 

0.62  X 10  ^ s 


The  characteristic  frequency  of  the  transient  is  (A/) 


16  kHz. 
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Appendix  B 

COMPARISON  OF  ENERGY  DENSITY  OF  COSMIC  PARTICLES 
AND  ENERGY  DENSITY  OF  AMBIENT  NOISE  IN  THE  OCEAN 

Ambient  noise  in  the  ocean  is  reported  in  decibel  units  referenced  to  the  intensity  (or 
energy  density)  of  a plane  wave  of  unit  amplitude  in  a 1-Hz  band.  Thus  if  is  the  energy 
density  actually  measured  and  ^re!  is  the  reference,  then  the  number  of  measured  decibels  is 

A = 10  log|o 

Hence 

f ^ = Erfi  10  (units:  ergjcm^  ■ Hz). 

For  example  at  25  kHz  the  lowest  measured  ambient  noise  is  reported  as  A = -82  dB  "re  1 dyne 
/cm  This  means  that 

= E,,!  10  2 = 6 3 X 10  -9 

Now  the  energy  density  of  a plane  wave  of  unit  amplitude  is 

- A 9)^  .4,44  , ,0  -,l 

pA  lx(1.5xl0^)^  cm^ 


from  which 


=4.44  X 10  ""  X 6.3  x 10  =2.8  x 20 


Using  the  same  procedure  we  construct  the  following  table. 


Freq.  (Hz) 

N (dB) 

(£^  (erg/cm  ^ • Hz 

1 

-11 

3.5  X 10 

10 

-40 

X 

o 

100 

-60 

4.4  X 10 

1,000 

-73 

2.2  X 10 

10,000 

-80 

4.4  X 10 

25,000 

-82 

2.8  X 10 

30,000 

-84 

1.76  X 10 

Let  us  assume  that  for  various  reasons  we  are  requied  to  use  a 10-kHz  band  for  underwater 
detection  of  these  particles.  The  average  noise  energy  density  in  a 1-Hz  band  between  10,000 
and  25,000  Hz  is  ~ 3.6  x 10  (erg/cm  ^ Hz).  Thus  in  a 10^-Hz  band  the  average  energy 
density  is 

3.6  X 10  X 10“*  =3.6  X 10  "'^erg/cm^ 

We  desire  to  compare  this  with  the  energy  density  of  the  cosmic-particle  flux.  To  do  this,  we 
let  /j,  be  the  particle  flux  density  per  steradian  and  E^  be  the  average  energy  per  particle. 
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Then  the  total  energy  flux  (intensity)  is  obtained  by  integrating  over  the  solid  angle  of  in- 
cidence; 


/^,  = J /jj  Ep  (HI  (units:  erg/t/n^  • s). 


We  next  assume  that  the  particle  flux  is  a plane  wave  of  light.  We  then  estimate  the  energy 
density  E to  be 


£■„  = 


Ip  S In  Ep 


(units:  erg/cm  ^ ), 


in  which  's  the  speed  of  light. 

As  an  estimate  we  take  the  case  of  muons  and  choose  /jj  =1.83  x 10  “^/cm’  ■ s • sr 
X cos  ^Oandf^,  = 2 xl0’eV=2  x 10^  xl.6  xl0“’^=3.2  x 10“^  erg/muon.  Over  the 
hemisphere  of  incidence 


1 .83  X 10  ^ muons  ( 1 f 


cm  • s • sr 


f cos^  0 277  sin  0 (/O  3.2  X 10  ^ erg/n 

•b 


= 1.83  X 10  3-2  X 10  -^  = 1.23  x 10  erg/cm  ^ s. 

The  energy  density  of  the  equivalent  plane  wave  of  "light"  is 

^1.23  X j_0-^erg_/cm^  ^ ^ ^ 10  "'5^. 

3 X 10 cm/s  cm-’ 

Thus,  if  we  listen  underwater  over  a 10-kHz  band,  the  energy  density  of  the  muon  flux  is 
about  the  same  as  the  energy  density  of  noise  in  the  ocean  over  the  band  of  10.000  Hz  to 
25,000  Hz.  On  the  other  hand,  if  we  listen  over  a 1-Hz  band,  the  noise  due  to  the  muous  (of 
energy  10^  eV)  is  about  10,000  times  higher  than  the  ambient  noise  due  to  all  other  causes. 
However  the  presence  of  a muon  signal  is  a random  event  in  space  and  time.  The  probability 
of  detection  is  still  to  be  calculated. 


Appendix  C 

TOTAL  NOISE  OF  MUONS 


r 


1 


Assume  that  muons  of  energy  2 GeV  penetrate  the  ocean  with  a total  (lux  of 
1,82  X 10  'Vcm^  ■ s • sr.  The  angular  distribution  of  particles  varies  as  cos^  where  H is  the 
angle  of  the  muon  track  with  the  vertical.  The  rate  of  absorption  is  2 MeV/cm,  so  that  the 
effective  range  of  the  track  is  Z.  = 10  m.  We  wish  to  calculate  the  total  noise  production  of  the 
muons  and  compare  it  with  the  noise  of  molecular  agitation  in  the  ocean.  To  do  this,  we  use 
the  formulas  of  model  la.  For  kL  large  the  total  acoustic  power  generated  by  vertically  incident 
particles  is 

« w (/„  S)^ 

^ 1 6 /> 

First  we  a,ssume  that  a = = (1/10)  m Then  we  must  replace  the  uniform  intensity  /„ 

of  vertically  incident  particles  by  an  effective  intensity  to  account  for  the  cos^  w dependency  of 
arrival.  As  a simple  approximation  we  lake 

IT  - 

'ilf  = 4"  / A,  f 277  sin  » M = ~^ 

Next  we  calculate  the  input  heat  power  IF,  to  the  water  caused  by  one  muon,  which  according 
to  the  model  is  to  be  based  on  the  heat  deposition  time  T,: 

^1  = ^ ^ ^ A.  ^ ^ ^ (units;  A • mfs). 

To  choose  7,  we  assume  first  that  the  noise  caused  by  muons  is  not  greater  than  the  measured 
ambient  noise  of  the  ocean  at  the  selected  frequency.  As  we  shall  see  later,  the  smallest  7,  = 
41  ns.  Thus 

X 2 GeV  X 1.6  X iQ  N • m/e  f 
' 3 ^ 0.41  ns 

= 1 .635  X 10  N • m/s. 


The  total  acoustic  power  "radiated"  by  one  muon  at  a frequency  of  25  kHz  (say)  is  therefore 

1 o>  )3  2 ( W,  ) 2 
= -L 


16  ,,  k-j 

= X 277  X 25  X 10-’  s X X 

10  m 16 

^ 11.1  > , 10 


1 


10^  N ■ s‘/m 


2/„,  4 


1 s^K' 


(4,18  X 10-’)^  m'* 


= 2.94  X 10 


-19 


N ■ m 


s muon 

We  are  given  that  there  are  1 .83  x 10  muons/s  • cm  A or  183  muons/s  m so  that  the  to- 
tal power  of  A muons  is 
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^^  V =2.94  X 10 


= 5.39  X 10 


J±_r^  ^ 183 


N ■ ni 


We  inlerprel  IV'y  as  the  acoustic  power  (watts)  flowing  through  a spatially  averaged  1 of 
ocean  surface  temporally  averaged  over  1 s.  We  will  go  one  step  farther  and  assume  the  acous- 
tic power  so  generated  belongs  to  a plane  wave  of  sound  arbitrarily  oriented,  the  energy  density 
of  which  is 


c = speed  of  sound  in  water. 


E x JO  -'^N  n,/s  ■ m^  ^ 3 ^ 

1 .5  X 10  ^ m/s 

= 3.59  X 10  erg/cm  ^ (averaged  over  1 s). 


-JO  N • m 


This  is  nearly  the  energy  density  of  ambient  noise  (per  herz)  measured  in  the  ocean  between 
15  and  25  kHz.  The  equality  of  muon  noise  and  ambient  noise  was  deliberately  made  to  occur 
by  choosing  T^  to  be  41  ns.  The  energy  density  is  calculated  on  a 1-s  averaging  basis  as  re- 
quired by  the  specified  incidence  rate. 

The  input  energy  £„  is  actually  taken  from  the  spectrum  of  measured  energy  versus  in- 
cidence rate.  It  is  an  arbitrary  selection.  Similarly  T^,  the  time  factor  in  the  rate  of  deposition 
of  heat,  is  taken  to  be  a plausible  estimate  in  the  absence  of  concrete  experimental  fact.  Previ- 
ously lor  various  reasons  we  adopted  T|  to  be  1 ns.  We  will  now  make  a calculation  based  on 
this  choice. 


For  T]  = 1 ns  the  heat  deposition  of  a single  muon  is 

11/  — 217  2 X GeV  X 1.6  x 10  N • m/eV 


(Ft  X 


= 0.67  N ■ m/s. 

The  acoustic  power  at  25  kHz  of  one  muon  is  therefore 


. 2,  X X -1 

10  m s 1< 


(1.4  X 10 


>6  10^  N • s/m^ 


X (6.7  X 10  -■) 


-uz  ■ m^ 


I s^  • K' 


s^  (4.18  X 10^)2  m'* 

= 4.95  X 10  N ■ m/s. 

We  multiply  this  by  the  incidence  rate  to  obtain  the  unit  areal  power  flux: 


=4.95  X 10 


-16  N • m 183  muons 


= 9.05  X 10  -'■♦Jl—Ql  X ( — i 


m ^ ■ s 
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This  is  then  divided  by  the  speed  of  sound  c to  obtain  the  space-and  time-averaged  energy 
density  of  an  equivalent  plane  wave: 


£ = 9.05  X 10  N— HI  X ( — ~ 


\ 1 / N ■ m 10 

-)  X ; =6.04 ; — X 

1.5  X 10^  m/s  m^  s 


= 6.04  X 10  erg/cm  ^ (averaged  over  1 s). 

This  is  larger  than  the  ambient  noise  in  a 1-Hz  band  by  a factor  of  6.04  x 10  “’V3.59  x 
10  = 1.7  X 10^.  It  is  seen  to  be  10  x smaller  than  the  noise  energy  density  in  a 10-kHz 

band. 

From  this  calculation  shows  that  the  prediction  of  noise  generation  from  muons  in  the 
ocean  is  highly  sensitive  to  the  correct  estimate  of  heat  deposition  rate.  Since  an  upper  limit  of 
this  noise  (at  25  kHz)  is  the  measured  ambient  noise,  one  can  find  for  any  choice  of  £„ 
with  the  understanding  that  it  is  the  smallest  possible  time  in  the  limit  that  all  the  ambient 
noise  is  due  to  muons  of  this  energy. 

A final  calculation  based  on  the  experimentally  more  plausible  value  of  Ti  ~ 1 ns  will 
serve  to  illustrate  the  wide  divergence  in  predictions.  Thus 

27r  2GeV  x 1.6  x 10  N • m/eV 

— X - 

‘3  \ fxs 

= 0.67  X 10  N ■ m/s 


Thus  the  "particle  power"  which  is  considered  in  this  model  to  be  the  analog  of  the  laser  power 
of  Westervelt-Larson  is  0.67  mW.  The  total  acoustic  power  radiated  by  one  muon  at  25  kHz  is 
then 


1 . 277  X 2.5  X lO''  (1.4  X 10 

10m  s 


X (0.67  X 10  -3)2 


1 1 V.  7,  4 Is^'K^ 

X — X N • s 7m  ^ X — - — - 

16  10^  (4.18  X 10^)^  m'^ 

= 4.94  X 10  “32  N . m/s. 

For  an  incidence  rate  of  183  muons/s  • m^  the  power  generated  over  1 m^  in  an  averaging 
time  of  1 s is 

=4.94  X 10  -32  . m/s  muon  x 183  muons 

= 9.04  X 10  -3®  watts  of  acoustic  power. 

Considering  again  each  square  meter  to  be  a unit  area  of  a plane  wave,  we  find  the  energy  den- 
sity to  be 

£ 9.04  X 10  -30  N ■ m/s  ■ m3  ^ N ■ m 

1.5  X lO^m/s  m3 

= 6.03  X 10  -32  erg/cm  3, 

The  ratio  of  muon  noise  to  ambient  noise  (in  a 1 Hz  band)  is  then 

_ 6.03  X 10  -32  ^ ^ 

^ambient  3.59  x 10 
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The  comparison  with  molecular  agitation  noise  is  different.  According  to  Mellen  (14]  the 
thermal  noise  in  the  ocean  is  empirically  given  by 

A,/,  - -175  + 20log|o  / 

in  units  of  dB  re  1 ^bar  of  equivalent  plane  wave  of  unit  amplitude,  with  / being  given  in 
hertz.  In  arithmetic  units  the  energy  density  of  thermal  noise  is  therefore 

£ =£  I / y 

^'■'  15.623  X 10«  I ■ 

Here  we  take  the  reference  energy  density  to  be  that  of  a plane  wave  of  unit  amplitude: 

£ (1)^  (dyne/cm _erg_ 

pc^  1 dyne_^  ^ 5 ^ jq5)2  cmi  cm^ 

cm'*  s^ 

Thus  the  thermal  energy  density  is 


E,„  =4.44 


X 10 


-11 


./ 


5.623  X 10* 


erg 
cm  ^ 


At  / = 25  kHz, 


E,„  =4.44 


X 10 


2.5  X 10'* 


5.623  X 10* 


= 8.78  X 10 


cm 


3 ■ 


The  ratio  of  acoustic  energy  density  caused  by  muons  to  the  noise  energy  density  of  molecular 
agitation  depends  on  the  bandwidth  of  listening  and  on  the  choice  of  deposition  time  T^.  Us- 
ing the  results  just  calculated,  we  can  construct  a ratio  (Table  Cl).  Thus,  to  detect  muon  noise 
above  thermal  noise,  the  bandwidth  of  the  receiver  must  not  be  appreciably  more  than  1 Hz. 


Table  Cl  — Ratio  of  25  kHz 


Bandwidth 

(kHz) 

T^  = 41  ps 

Tj  = 1 ns 

1 

10 

3.59  X 10-'9/8.78  X 10-^0=*=  4 

3.59  X 10-'9/8.78  X IQ-**  ^ 4 X lO"^ 

6.04  X 10-'®/8.78  X 10*20  7 x 10^ 

6.04  X 10-'6/8.78  X lO*'*  « 0.7 

REMARKS  I 

The  procedure  used  in  the  above  calculation  is  to  calculate  the  acoustic  radiation  from 
one  muon  and  then  use  the  statistics  of  muon  incidence  to  calculate  the  power  per  unit  of  time 
per  unit  of  irea.  The  statistical  basis  of  the  muon  count  must  not  not  be  overlooked.  Most 
likely,  muon  incidence  corresponds  to  shot  noise  and  has  therefore  a Poisson-type  probability 
distribution  in  time.  Over  a long  record  (say  many  seconds,  hours,  days,  etc.)  the  average  of 
the  distribution  is  taken,  and  this  is  the  number  used  in  the  above  calculation.  Other  probabil- 
ity moments  tan  be  of  use  in  developing  the  statistics  of  acoustic  radiation  developed  by  the 
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muons.  For  example,  by  quoting  the  variance  of  the  random  incidence,  we  can  find  the  vari- 
ance of  the  power  crossing  a square  meter  of  the  ocean  surface;  hence  we  can  find  the  vari- 
ance of  the  acoustic  energy  density.  This  would  be  a valuable  statistic  in  comparing  the  acous- 
tic radiation  developed  by  muons  to  the  noise  of  molecular  agitation  in  the  ocean. 

We  have  made  our  calculation  on  the  basis  of  a component  (25  kHz)  of  the  power  spec- 
trum of  the  radiated  acoustic  noise,  choosing  a single  muon  particle  as  our  basic  model.  When 
we  deem  iV  muons  incident  in  1 s (as  averaged  over  many  seconds  of  record),  we  have  as- 
sumed incoherent  summation  by  taking  the  total  power  to  be  yv  times  the  power  of  one  muon. 
If  the  summation  is  coherent,  there  will  be  a directional  gain  in  radiation,  as  if  the  volume  of 
incidence  were  a volume  array  of  sources.  The  energy  density  would  then  be  higher  in  specific 
directions. 

We  can  also  treat  the  problem  in  the  time  domain  rather  than  in  the  frequency  domain. 
Thus  we  can  calculate  the  transient  power  tV(t)  per  muon,  find  the  average  power 


= r dt/M 

*0 

for  the  duration  It  of  the  transient,  and  then  multiply  by  the  statistically  determined  time  and 
space  rate  of  incidence  S to  find  the  total  power  for  all  frequencies  per  unit  of  time  and  area. 

Finally  we  could  proceed  to  calculate  a long  time  record  of  power  IF(r),  taken  over  many 
seconds,  and  subject  it  to  a probability  analysis.  This  can  be  done  by  dividing  the  record  into 
M records  and  making  histograms  of  the  average,  the  variance,  etc.  among  the  components  of 
the  set.  An  autocorrelation  could  then  lead  directly  to  the  power  spectrum. 

REMARKS  II 

The  total  noise  energy  density  generated  by  the  flux  of  muons  is  an  average  over  all 
space.  The  sharp  directionality  of  vertically  incident  particles  has  been  canceled  in  the  averag- 
ing process.  Although  illustrative,  the  calculation  of  acoustic  energy  densities  does  not  lend  it- 
self to  experimental  check  directly.  One  must  measure  acoustic  pressure  and  use  it  to  calculate 
energy  density  according  to  some  simplifying  assumption,  say  plane  waves.  The  significant 
directional  power  quantity  is  the  acoustic  intensity.  From  Model  la  in  the  main  text  this  inten- 
sity is  seen  to  depend  on  the  angle  ft  as 

1 

1 + (kL)^  cos^  ft 

in  which  ft  is  measured  relative  to  the  vertical.  Actually  muons  arrive  at  all  angles  "£  (relative 
to  the  vertical).  The  sound  radiated  from  a track  at  angle  ( in  direction  of  angle  i from  £ then 
depends  on 

cos^  £ 

1 + ikD^cos^  {£  -f  ^)  ■ 

Thus  the  greatest  contributions  to  muon  noise  occur  in  the  horizontal,  from  tracks  near  the 
vertical  (f  — 0).  A rough  approximation  to  noise  generation  is  to  average  cos^  £ in  the 
numerator  over  a solid  angle  of  hemisphere  and  to  neglect  £ in  the  denominator. 
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REMARKS  III:  DIRECTIONALITY 

The  number  of  particles  incident  on  a unit  area  is  proportional  to  cos^  W,  where  H is  the 
angle  of  the  track  of  the  particle  relative  to  the  vertical.  Since  the  radiated  sound  from  a track 
is  sharp  in  a direction  normal  to  the  track,  we  assume  that  the  directionality  of  all  the  tracks 
through  a unit  area  is  also  proportional  to  cos^  W,  where  H is  now  taken  relative  to  the  horizon- 
tal. The  3-dB  points  of  the  miensiiy  pattern  of  acoustic  radiation  occurs  at  H values  such  that 
COS'  ft  =1/2,  or  cos  H =0.707,  or  « = ± 45°  . Hence  the  3-dB  beamwidth  of  muon  noise  is 
90°  centered  on  the  horizontal.  We  compare  this  directionality  with  that  of  ambient  noise  and 
of  directional  noise.  Ambient  noise  at  low  frequencies  (say  less  than  10  kHz)  normally  arrives 
in  a wedge  ±20°  centered  on  the  horizontal  [Ref.  22,  p.  282).  High-frequency  ambient  noise 
originating  in  surface  wave  action  arrives  at  a point  in  a vertical  cone  whose  vertex  angle 
depends  on  depth  of  the  listening  sensor.  In  the  range  of  interest  (25  kHz  and  above)  high- 
frequency  muon  noise  can  be  differentiated  from  ambient  noise  by  the  difference  of  angle  of 
arrival:  muon  noise  arrives  horizontally,  and  surface  noise  arrives  vertically.  In  contrast  ther- 
mal or  molecular-agitation  noise  dominates  at  high  frequency  and  is  omnidirectional. 

REMARK  IV:  RELATION  BETWEEN  ENERGY  DENSITY  AND  INTENSITY  14,  p,  5771 

We  imagine  the  total  acoustic  effect  of  cosmic  particles  in  the  ocean  to  be  analogous  to 
the  assemblage  of  randomly  distributed  standing  plane  waves  in  a room,  each  wave  traveling  in 
some  direction  <t>,  ft.  Let  each  plane  wave  have  amplitude  A (r,  <t>,  w).  The  power  incident  on 
a unit  area  normal  to  <A,  W is  the  intenstiy  /: 

. 2tt  TT.'2 

/(r ) = — r d<b  r \A{r,  cos  (t  sin  H dH. 

f>c  i i 

Similarly  the  energy  density  w (r)  is, 

27r  T 

fV(r)  = — j J d<t>  J I /I  (r,  <i>,  W )|  ^ sin  t)  cM. 

f>C^  0 *0 

At  high  enough  frequency  we  can  take  the  sound  field  in  the  room  to  be  isotropic,  meaning 
that  A ir,  <t>,  H)  is  independent  of  r,  <fr,  «.  This  is  the  frequency  range  of  geometrical  acoustics. 
Upon  integrating  in  this  region,  one  obtains 

, A^  TT  Att  A ^ 

I = , W — . 

pc  pc  2 

Hence  the  relation  of  / to  w is. 


In  the  above  numerical  calculation  of  muon  noise  we  have  omitted  the  factor  4 and  written  w 
= He.  This  omission  was  done  because  an  integration  over  the  hemisphere  of  particle  in- 
cidence had  already  converted  angular  incidence  to  vetical  incidence. 

STATISTICS  OF  NOISE  IN  THE  OCEAN 

The  power  spectrum  level  of  noise  in  the  ocean  has  been  complied  by  Wenz  [13]  and  is 
reproduced  as  Fig.  Cl  for  convenience.  The  ordinate  of  this  chart  represents  the  noise  level 
(AT),  defined  as 
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10^  103 

FREQUENCY  CPS 


I'ig.  C'l  - Wcnz  chart  of  noise  in  the  ocean 
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NL  — \0  log  10  P f = frequency  (//;) 

Pi,  (/) 

in  which  the  reference  is  the  quantity  of  a plane  wave  of  unit  amplitude  at  the  frequency  j. 
All  curves  shown  (except  the  one  labeled  molecular  agitation)  represent  measured  values  of 
ambient  noise.  The  curve  for  molecular  agitation  is  a straight  line  extending  to  the  right  with 
the  slope  shown.  All  quantities  are  based  on  noise  power  in  a 1-Hz  band.  The  lower  limit  of 
ambient  noise  in  the  ocean  at  25  kHz  is  -8  -74  = -82  dB  re  1 dyne/cm  ^ (spectrum  level), 
and  the  noise  of  molecular  agitation  is  -13  -74  = -87  dB  re  1 dyne/ cm  By  use  of  the 
Mellen  formula  [14]  one  can  extend  this  latter  curve  upward  as  desired. 

Signal  Processing  of  Muon  Noise  When  the  Time  of  Arrival  is  Random 

We  form  a conceptual  analog:  let  the  unit  area  of  the  ocean  surface  be  likened  to  an 
anode  in  a vacuum  tube,  and  let  the  cosmic  particles  incident  on  it  be  likened  to  electrons  ar- 
riving randomly.  Also  let  each  particle  generate  the  same  heat  source  h(x,  i)  in  the  water  and 
hence  the  same  acoustic  pressure  pit,  i).  Assume  all  effects  of  a unit  area  of  incidence  are  ad- 
ditive, taking  r the  same  for  all.  The  total  effect  at  time  t due  to  all  the  particles  of  a small 
enough  unit  area  is 

oo 

P(r.  ')  = Z p(t,  I - ). 

k = 

Since  the  individual  occurs  at  random  time  the  total  Pis  a random  function  of  time. 

Choose  an  interval  T and  let  exactly  K cosmic  particles  arrive  on  the  unit  area  in  this  in- 
terval at  random  times  Then  the  contribution  of  these  K particles  to  the  total  time  history 
of  acoustic  pressure  is  called  Py^-,  where 

*. 

Pk  = L /’<'■• ' “ 'a  >• 

k -I 

If  r is  dropped  for  convenience,  the  autocorrelation  <)>  (r)  of  the  total  acoustic  pressure  will  be 


i//  (r  ) = Pit)  p(r  + t)  = ^ piK)  Pi^  it)  p,^  (l  + r ), 

k -0 

in  which  piK)  is  the  probability  that  exactly  K muons  arrive  in  the  interval  0,  T.  If  one 
chooses  piK)  to  be  a Poisson  distribution  corresponding  to  shot  (or  impulse)  noise,  then  one 
can  show  that 

oo 

ill  (t)  = a J"  pit)  pit  + t)(Ii  + pit)^. 

— oo 

in  which  A is  the  average  number  of  muons  arriving  per  second  and  Pit)  is  the  average  value 
of  Pit),  namely, 

oe 

pit)  = A J"  pit)  dt. 

— CO 

The  ensemble  averaged  power  spectrum  iv(/)  can  now  be  calculated  from  the  standard  relation 

oo 

w(/)  =2  J*  ill  ir)  cos  277  / t dr 

— oo 

= 2A|,v(/)|  2 -f-  2 PTtVfiif) 


I 


i 


i 


J 
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in  which 

oo 

s{/)  = / p(i)  e dt. 

— OO 

Thus,  to  find  the  power  spectrum,  one  must  Fourier  transform  the  impulse  acoustic  pressure  of 
one  muon  to  obtain  the  amplitude  spectrum  s(J).  Then  one  must  take  2A^ times  the  absolute 
value  squared  of  s(f)  and  finally  add  to  it  twice  the  square  of  the  mean  value  of  the  total  pres- 
sure P(t)  as  a "DC"  component. 


When  w(f)  is  available  over  the  frequency  range  of  interest,  we  can  estimate  the  statisti- 
cal behavior  of  PU).  A convenient  reference  is  S.  O.  Rice  ("Mathematical  Analysis  of  Random 
Noise,"  in  Selerected  Papers  on  Noise,  N.  Wax,  editor,  Dover  Publications,  1954).  As  a special 
case  we  list  the  applicable  formulas  for  narrowband  noise,  assuming  the  total  muon  noise  P(t) 
will  be  processed  through  a narrowband  filter  to  obtain  spectral  information  at  (say)  25  kHz. 


Thus  assume  the  autocorrelation  is  of  the  for  m i1<(t)  = A cos  fir-,  then,  if  P(!)  and 
BP/di  are  both  Gaussian,  the  expected  number  of  zero  crossings  per  second  is 


<M0) 


<i(0) 


1/2 


= 2 


J./^  >v(/)  df 


f w{f')  df 


1/2 


where 


</»''(0)  = — 4Tr^  r.T^»»'(/')  df. 

The  probability  that  the  time  interval  between  successive  zeros  lies  between  t and  t + dr  is 
approximately 

dr a 

2 tl  -1-  (t  — T|  ) ^)  ^^^ 

where 


a 


ifb 

h - fa 


and  T| 


in  which  — /^  is  the  passband. 


1 

h + fa  ' 


Writing  P(t)  — cos  — Pj  sin  w,,/,  where  (o„,  = 27r/,„  is  the  midfrequency  of  the 
filter  and  P,.  and  Pj  are  quadrature  components,  we  take  P^.  and  P^  as  normally  distributed  with 
the  same  standard  deviation.  The  statistics  of  the  pair  P^.  and  P,,  and  of  the  envelope  P(t)  = 
-I-  Pj  are  the  following: 

• the  probability  that  the  point  (P,.,  Pj)  lies  in  the  elementary  rectangle  dP^  dP.^  is 

dPc  dP.  P}  + Pi 

2Tr.^(0)  ' 
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• the  number  of  maxima  of  the  envelope  R(t)  expected  in  1 s is  0.6  (/^  — . 


• the  probability  that  a maximum  selected  at  random  lies  between  R and  i 
the  envelope  R > 2.5  Vi/TToT,  is 

cfR  R 


1.13  02  - i)  e 


the  expected  number  of  maxima  per  second  is 


y = 


] df 


J /^wO)  df 


1/2 


L 


fah 

? + dR,  when 


i 
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Appendix  D 

MAGNITUDE  OF  CONSTANTS  USED  IN  THE 
NUMERICAL  CALCULATIONS 


^ = 1.4  X 10  ^ °K  coefficient  of  thermal  expansion  of  water 

C =4  18  X 10 3 Qf — , specific  heat  of  water  at  constant  pressure 

f A>  °A  ^ec2  °K 

r^L  = 1.6  X 10  “'^joules  =1.6  x 10  erg,  energy  ^ 

£)  = 1 43  X 10  = 1.43  X 10  heat  diffusion  constant  of  water 

sec  sec 

C'  = 1.5  X 10  ^ speed  of  sound  in  water 
sec 

■ A/  dynes 

1 atmosphere  = 10^  — - = 10° — 

nr  cm'^ 

reference  pressure  on  dB  scale: 

(1)  dB  re  1 dyne/cm  ^ 

(2)  dB  re  0.0002  dyne/cm  subtract  74  dB  to  convert  to  dB  re  1 dyne/cm  ’ 

(3)  dB  re  1 micropascal,  subtract  100  dB  to  convert  to  dB  re  1 dyne/cm  ^ 

1 calorie  =4.18  joules  =4.18  x 10^  ergs  =2.6  x 10'‘^eV 
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